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Tlie  n  most  vital  links  of  a  flow  network  are  defined  as  those 
n  arcs  wliose  simultaneous  removal  from  the  network  causes  the  greatest 
decrease  in  the  throughput  capability  of  the  remaining  system  between 
a  specified  node  pair.   These  n  arcs  are  shown  to  be  tlie  n  largest 
capacity  arcs  in  some  cut.   A  solution  procedure  is  developed  which 
involves  sequentially  modifying  the  network  in  such  a  manner  as  to 
make  this  cut  minimal.   An  algorithm  with  computational  results  is 
presented. 


CHAPTER   1 


INTRODUCTION 


The  research  presented  in  this  dissertation  is  concerned  with 
the  problem  of  optimally  removing  components  of  a  flow  network.   There 
are  many  problems  embedded  in  this  part  of  network  theory.   This  is 
due  to  the  multitude  of  considerations  which  are  involved  in  stating 
and  solving  a  network  problem.   For  instance,  one  must  note  the  type 
of  network  involved;  which  and  how  many  network  el  ei'ients  are  candidates 
for  removal;  and  what  is  the  measure  of  optimality  (i.e.,  what  was  the 
purpose  of  posing  the  problem). 

The  problem  of  removing  arcs  from  a  network  so  as  to  stop  all 
flow  from  a  set  of  sources,  S,  to  a  set  of  sinks,  T,  at  a  minimal  cost 
(where  each  arc  requires  an  assigned  cost  for  removal)  is  referred  to 
in  the  literature  as  the  Attacker's  Problem.   The  related  problem  of 
reinforcing  a  network  so  as  to  most  increase  the  cost  of  the  optimal 
solution  to  the  Attacker's  Problem  is  referred  to  as  the  Defender's 
Probleiri.   In  this  latter  problem  there  is  an  assigned  cost  which  is 
required  to  improve  any  link  and  a  budget  consideration  which  limits 
the  "improvements"  that  can  be  made.   Thus,  there  are  several  sets  of 
improvenonts  which  the  defender  could  make  and  remain  within  his  budget. 
Any  such  set  of  improvements  is  called  cost  feasible  and  each  of  these, 
if  implemented,  would  lead  to  a  different  reinforced  network. 


The  Defender's  Problem  is  to  find  that  cost-feasible  set  of  improvements 
which  gives  a  reinforced  network  with  the  property  that    the  optimal 
solution  to  the  Attacker's  Problem  used  against  this  reinforced  network 
is  at  least  as  large  as  the  optimal  attacker's  solution  against  any  of 
the  other  feasible  reinforced  networks.   Jarvis  [19]  solves  both  the 
Attacker's  and  Defender's  Problems  for  the  single-comjnodity  case  by 
using  the  Generalized  Lagrange  Multiplier  Method  of  Everett  [10]  and 
Nunn  [24]  ,  and  the  transshipment  model  of  either  Busacker  and  Gowen  [6] 
or  Ford  and  Fulkerson  [11].   In  the  same  paper,  Jarvis  solves  the  multi- 
commodity  Attacker's  Problem  and  formulates,  but  does  not  solve,  the 
multi-commodity  Defender's  Problem.   Ratliff  [26]  solves  the  multi- 
commodity  Defender's  Problem  for  a  network  by  solving  a  series  of  multi- 
commodity  minimal  disconnecting  set  problems  and  a  series  of  set -covering 
problems.   Ratliff  uses  a  method  developed  by  Bellmore,  Greenberg,  and 
Jarvis  [3]  to  solve  the  former  series  of  problems  and  a  method  of 
Bellmore  and  Ratliff  [4]  to  solve  the  latter  series  of  problems. 
Greenberg  [15]  defines  a  more  general  Attacker's  Problem  than  that  of 
Jarvis.   His  method  of  solution  for  the  single-commodity  case  involves 
the  Generalized  Lagrange  Multiplier  Method  and  the  maximum  flow  algo- 
rithm [11].   He  solves  the  multi-commodity  Attacker's  Problem  by  modi- 
fying a  0-1  integer  progranmiing  method  described  by  Geoffrion  [14]  and 
Balas  [1]. 

It  can  be  noted  that  by  setting  the  cost  of  destroying  an  arc 
equal  to  its  capacity, the  Attacker's  Problem  becomes  synonymous  with 
the  minimal  disconnecting  set  problem.   Bennington  [5]  gives  a  metliod 
to  solve  the  latter  problem  for  an  undirected  network.   He  also  shows 
that  the  multi-commodity  minimal  disconnecting  set  problem  can  be  posed 


as  a  node  isolation  problem.   A  discussion  of  the  isolation  problem  and 
a  method  to  solve  it  are  given  by  Bellmore,  Bennington,  and  Lubore  [2]. 

The  notion  of  assigning  values  to  different  elements  of  the  net- 
work has  also  been  employed  by  various  people.   Jlobbs  [18]  assigns  num- 
bers to  the  links  and  nodes  of  a  communication  network  to  reflect  their 
relative  importance.   Tlie  determination  of  these  "im))ortance  numbers" 
which  are  based  on  the  configuration  of  the  network  and  the  assigned 
value  of  each  commander  —  subordinate  pair ,  is  equal  to  a  sum  of  products. 
The  sum  is  taken  over  the  possible  commander — subordinate  pairs,  where 
each  element  of  the  sum  is  a  product  of  the  value  of  the  particular 
commander — subordinate  pair  and  the  percent  of  chains  connecting  the 
commander  and  subordinate  which  contain  the  link  in  question.   Krokowski 
[20]  also  suggests  that  importance  values  be  assigned  to  links  and  nodes, 
but  he  proposes  that  a  good  measure  of  the  value  of  a  network  component 
(in  a  war  situation)  would  be  the  replacement  cost  of  the  element.   He 
recommends  a  linear  programming  approach  but  does  not  develop  the  solu- 
tion beyond  the  conceptual  stage. 

The   effect  of  changing  the  length  of  a  link  on  some  or  all  of 
the  shortest  paths  between  nodes  in  a  network  is  addressed  by  Haider 
[17].   In  fact,  he  locates  the  shortest  path  distances  that  will  be 
affected  by  increasing  the  length  of  an  arc  for  a  symmetric  distance 
network  (a  network  where  d(a,b)  =  d(b,a)  for  all  nodes  a,  b  in  the 
set  N). 

It  should  also  be  noted  that  the  general  subject  of  network 
vulnerability  and  synthesis  can  be  viewed  strictly  as  a  graph  theory 


p rob] em.   For  instance,  Frank  and  Frisch  in  [13]  present  analysis  meth- 
ods for  problems  similar  to,  and  yet  different  from,  those  previously 
mentioned.   In  i^articular,  they  consider  an  Attacker-like  Problem  in 
whicli  the  objective  is  to  determine  tlie  minimum  number  of  nodes,  arcs, 
or  arcs  and  nodes  that  must  be  removed  to:   (1)  stop  all  possible  flow 
between  designated  node  sets,  and  (2)  partition  the  network  into  two 
disjoint  sets  (i.e. ,  source  and  sink  nodes  are  not  specified  and  break- 
ing the  communication  between  any  pair  of  nodes  is  acceptable).   They 
also  consider  the  problem  of  constructing  a  network,  subject  to  cost 
consti-aints ,  for  which  the  number  of  elements  that  must  be  I'emoved  to 
partition  the  network  into  two  components  is  maximal. 

Problem  Definition 

The  problem  to  be  addressed  in  this  paper  is  to  determine  those 
n  arcs  (or  links)  wliose  simultaneous  removal  from  the  network  causes  the 
greatest  decrease  in  the  maximal  possible  flow  between  some  specified 
node  pair.   This  problem  will  be  referred  to  as  the  n  Most  Vital  Links 
Problem  and  the  best  set  of  n  arcs  to  remove  will  be  called  the  n  most 
vital  links.   Note  that  if  removing  n  arcs  stops  all  the  flow  between 
a  source  and  sink,  then  this  problem  becomes  the  Attacker's  Problem 
where  the  cost  to  destroy  any  arc  is  the  same. 

The  following  definitions  and  notation  will  be  used  throughout 
this  dissertation. 

A  network  will  be  denoted  by  G(N,A),  where  X  =  (x  ,x  ,.,.}  is 
a  finite  collection  of  elements  x. ,  each  representing  a  node  in  the 


network,  and  A  is  a  subset  of  the  j)airs  (x,y)  of  elements  from  N  repre- 
senting arcs  in  the  network.  (Another  characterization  of  the  arcs  which 
will  be  einiiloyed  is   a  eA.)  The  number  of  nodes  in  N  will  be  denoted  by 
[N|  and  the  number  of  arcs  in  A  by  m.   If  the  arc  (x,y)  is  an  ordered 
pair,  then  it  is  called  a  directed  arc;  otlierwise,  it  is  an  undirected 
arc.   For  the  results  presented  in  this  dissertation,  arcs  may  be  con- 
sidered as  either  directed  or  undirected. 

Given  the  network  G(N,A) ,  if  a  non-negative  number  denoted  by 
either  c(x,y)  or  c(a)  is  assigned  for  each  arc  in  A,  then  this  value 
will  be  defined  as  the  capacity  of  the  arc  (x,y),  as  it  can  be  t)iought 
of  as  the  maximal  amount  of  some  commodity  that  can  arrive  at  node  y 
from  node  x  pei-  unit  time. 

One  of  the  fundamental  notions  in  network  theory  is  the  concept 
of  the  maximal  flow  between  a  designated  pair  of  nodes  (s  and  t).  Given 
the  network  G(N,A)  and  a  capacity  c(x,y)  for  each  (x,y)  e  A,    a  question 
which  arises  is:   what  is  the  maximum  amount  of  steady  state  flow  that 
can  travel  from  node  s  to  node  t  over  the  network  G(N,A) ,  subject  to 
the  condition  that  the  flow  over  any  arc  is  no  larger  than  the  capacity 
of  the  arc?   Ford  and  Fulkerson  [11]  formulate  this  problem  as  a  linear 
programming  problem  with  a  special  structure.   This  structure  allows 
the  use  of  very  efficient  methods  (the  maximum  flow  algorithm)  for 
solving  the  problem.   Let   f (x,y)  be  a  variable  which  corresponds  to 
the  flow  over  arc  (x,y)  and  v  the  net  flow  out  of  node  s,  or  equiva- 

lently,  the  net  flow  into  node  t  (i.e.,  v  =   E  f(s,y)  -   Z  f(y,s) 

yeN         yeN 


foi'   all   pairs  which  are  arcs  in  the  network).   A  statement  of  the 
maximal  flow  problem  is: 


Max  V 

such  that 


if   X  =  s 


Z  f(x,y)  -  E  f  (y.x)  =  (  0    if  x  ?^  s,t 

yeN         yeN 

V -V    if  X  =  t 


f(x,y)  ^  c(x,y)   for  all  (x,y)  e  A 
f (x,y)  ^  0       for  all  (x,y)  e  A 

Finally,  two  related  sets  of  arcs,  disconnecting  sets  and  cuts, 
will  be  discussed.  A  cut,  (X,X),  separating  a  designated  pair  of  nodes 
s  and  t  is  a  set  of  arcs  defined  by: 

(X,X)  =  {(x,y)  I  (x,y)  eA,  xeX,  and  y  e  X} 

where  the  sets   X  and  X  are  any  partition  of  the  set  of  nodes,  N,  such 

that   s  ex,  t  eX,  X  f^,  X  =  (j) ,  and   XUX  =  N.  (A  cut,  on  occasion,  will  also 

be  denoted  by  A*.)  A  disconnecting  set  separating  a  designated  pair  of 

nodes  s  and  t  is  a  set  of  arcs  which  partitions  the  nodes  of  N  into 

p  sets:  X,  ,X„ , .  .  .  ,X   such  that  if  s  £  X  ,  then  t  ,6  X.  for  i  -    1,2,...  ,p; 
1   2      p  1  1 

X.  nx.  =  h      for  i  ?i'  i  :  and   U  X  =  N.   Thus,  a  cut  is  a  special  type 
1   J  i=l   i 

of  disconnecting  set,  i.e. ,  a  disconnecting  set  for  which  p  =  2.   The 
capacity  or  value  of  a  cut  and  a  disconnecting  set  will  be  taken  as  the 
sum  of  the  capacities  of  the  arcs  in  the  cut  and  disconnecting  set, 
respectively. 


Letting  a  cut  be  defined  as  minimal  if  its  value  is  at  least 
as  small  as  the  value  of  every  cut  of  G(N,A) ,  Ford  and  Fulkerson  [ ll] 
prove  that  the  maximal  flow  value  of  G(N,A)  is  equal  to  the  minimal 
cut  value  of  the  same  network. 

Using  the  notation  and  results  just  discussed,  the  n  most  vital 
links  problem  may  now  be  restated  as: 

Let  E.  for  i  -    1,2,.., ,p  be  all  subsets  of  A  that  contain 
exactly  n  elements  (the  cardinality  of  the  set  E.  will  be 
denoted  by  |e.|),  and  let  C  (N,A)  be  t]ie  capacity  of  a  min- 
imum cut  in  G(N,A).   Find  the  arc  set  E,  such  that 

k 

C  (N,A-E,  )  S  C  (N,A-E.)  for   j  =  1,2,...,  p. 
o       k     o      J 

Applications 

Some  applications  of  the  most  vital  links  problem  include: 

(1)  A  conflict  situation  where  there  is  a  logistics  or  communi- 
cations network,  a  defender  or  user  of  this  system  and  an  interdictor 
of  the  system.   The  user  wishes  to  know  which  arcs  are  most  vital  to 
him  so  that  he  can  reinforce  tliese  arcs  against  attack,  while  the 
attacker,  naturally,  wants  to  destroy  those  arcs  which  most  decrease 
the  efficiencj'  of  the  system. 

(2)  A  traffic  control  situation  where,  during  evening  rush 
hour,  most  motorists  go  from  the  city  (source)  to  the  suburbs  (sinks). 
A  question  which  arises  is:   If  n  road  segments  were  to  become  blocked 
what  effect  would  this  have  on  the  rush  hour  flow  and  the  blockage  which 


n  road  segments  would  most  decrease  flow?  One  use  for  this  information 
arises  in  the  scheduling  of  work  crews  to  the  transportation  network. 
For  example,  given  there  are  n  different  repair  crews,  simultaneously 
assigning  them  to  the  n  most  vital  links  would  not  be  advisable. 

Previous  Contributions 

The  most  vital  links  problem  and  other  very  similar  problems 
previously  have  been  addressed.   These  findings  are  mentioned  in  this 
section  with  the  most  significant  results  being  outlined  in  more 
detail. 

V/ollmer  [29]  has  devised  an  algorithm  for  finding  the  single 
most  vital  link  in  a  network.   Also,  by  solving  shortest  path  problems 
[11],  for  a  dual  graph  [7],  in  [27]  he  has  been  able  to  solve  for  the 
n  simultaneous  most  vital  links  for  a  planar  undirected  netv/ork. 
Durbin  [9]  uses  Wollmer' s  single  most  vital  link  approach  to  target 
strikes  for  an  attacker  of  a  given  network.   Each  arc  is  given  a  repair 
time  (in  time  periods)  and  n  strikes  are  made  for  each  time  period  on 
that  part  of  the  network  which  is  currently  operative.   Ondrasek  and 
Wollmer  [25]  address  a  problem,  similar  to,  but  more  general  than  that 
of  Durbin  as  they  add  tlie  consideiation  of  repair  time  and  cost  of  an 
arc  in  determining  whicli  arcs  to  attack. 


Tlie  single  most  vital  link  problem  allows  solution  methods 
which  cannot  be  efficiently  applied  to  the  !iiore  general  problem.   Foi* 
this  reasoii  c>nly  results  concerning  the  more  general  problem  are  pre- 
sented here.   Chapter  2  will  be  solely  devoted  to  the  single  most 
vital  link  problem. 


The  most  significant  contributions  in  this  area  are  due  to 
Wollmer  [27,28],  and  McMasters  and  Mustin  [22].   Their  results  can  be 
sunuaarized  as  follows: 

Results  due  to  Wollmer 

Letting  a  network  G(N,A)  be  defined  as  planar  if  it  can  be 

drawn  on  a  plane  such  that  no  two  arcs  intercept  except  at  nodes,  and 

as  source-sink  planar,  if  the  network  generated  by  adding  to  G(N,A)  an 

arc  from  the  source  to  sink  is  planar,  the  problem  proposed  by  Wollmer 

can  be  stated  as: 

Consider  a  source-sink  undirected  planar  network  in  which 
each  arc  is  subject  to,  at  most,  one  breakdown:   a  breakdown 
resulting  in  a  reduction  in  the  capacity  of  that  arc  by  a 
known  quantity.   It  is  desired  to  find  the  greatest  reduction 
in  maximum  flow  possible  from,  at  most,  n  breakdowns  and  where 
these  breakdowns  must  occur,  to  give  this  reduction. 

If  a  given  network,  the  primal  network,  is  source-sink  planar 

there  exists  a  related  dual  network.   Wollmer' s  solution  technique  uses 

these  two  networks  and  a  relationship  developed  by  Ford  and  Fulkerson 

in  [12];  i.e.,  the  result  is: 

A  one  to  one  correspondence  exists  between  the  cut  sets 
of  the  original  or  primal  network  and  the  loopless  path  through 
the  dual  network  (in  the  dual  network  the  arcs  have  lengths 
rather  than  capacities).   The  problem  of  finding  the  minimum 
cut  set  in  the  primal  network  is  equivalent  to  finding  the 
shortest  path  through  the  dual. 

Thus,  finding  the  maximal  flow  value  of  a  source-sink  planar 
network  is  equivalent  to  solving  a  shortest  path  problem  through  the 
dual.   Also,  the  problem  proposed  by  Wollmer  can  be  stated  in  terms 
of  the  dual  network  as: 


10 


Let  eacli  arc  of  the  dual  network  have  two  lengths  assigned 
to  it,  i.e. ,  the  original  and  the  reduced  length  (where  in  the 
primal  network  these  correspond  to  the  original  and  reduced 
capacities,  res])ectively) .   Assuming  that  the  lengths  of  at  most 
n  arcs  can  be  reduced,  the  problem  is:   the  length  reduction 
of  which  n  arcs  will  lead  to  the  shortest  possible  path  through 
the  network? 

Wollmer,  using  the  fact  that  shortest  path  problems  are  in  gen- 
eral easier  to  solve  than  other  network  problems,  developed  an  algo- 

th 
rithm  which  solves  n  +  1  network  problems.   The  first  or  0   problem  is 

a  shortest  path  problem  for  the  dual  network.   The  other  n  problems, 

for  i  =  l,...,n,  solve   the  following:   Given  that  the  shortest  path 

problem  in  which  at  most  i  -  1  arcs  may  have  their  length  decreased  has 

been  solved,  what  is  the  optimal  solution  to  a  shortest  path  problem 

v/hich  allows  at  most  i  arcs  to  have  a  reduced  length?   Wollmer  solves 

the  i    Eubproblem  by  finding  the  shortest  i-arc  paths  between  the 

source  and  all  the  other  nodes  of  the  network,  where  an  i-arc  path  to 

node  a  is  anj'  path  from  the  source  to  node  a  along  which  at  most  i  arcs 

have  their  length  reduced.   The  algorithm  uses  information  gained  by 

solving  the  i  -1    subproblem  to  solve  the  i    subproblem,  i.e.  ,  all 

the  shortest  i  -  1  arc  path  lengths  are  used  in  finding  the  set  of  the 

shortest  i-arc  paths. 

Results  due  to  McMasters  and  Mustin 

The  interdiction  problem  posed  by  McMasters  and  Mustin  in  [22] 
is  not  exactly  the  most  vital  linlis  problem  as  defined  in  this  paper. 
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However,  the  problem  which  they  solve  is  very  closely  related  to  the 
stated  problem  and  is  a  natural  parallel  problem  or  extension  to  the 
most  vital  links  problem.   Thus,  the  problem  statement  and  a  brief 
description  of  tJie  solution  method  aregiven  in  this  section. 
The  problem  in  the  notation  of  its  originators  is: 

Consider  an  undirected  planar  network,  [N;A].   Each  arc 
(i,j)  in  the  network  is  assumed  to  have  a  maximum  flow  capac- 
ity, u.  .  ^0,  and  a  minimum  flow  capacity,  1  ^.  0  ,    and  it  is 
assumed  that  at  least  one  arc  has  1  . >  0 ,   After  the  inter- 

diction  has  taken  place,  the  actual  capacity  of  an  arc,  m 

ij 

will  satisfy   0^1    ^  m   ^  u   . 
ij     iJ     ij 

Assuming  tliat  the  interdictor  incurs  a  cost,  C  .,  per 
unit  of  capacity  destroyed,  then  this  total  cost  for  reduc- 
ing an  arc's  capacity  from  u. .  to  m  .  is  C  .  [u    -  m   ]. 

ij      iJ      iJ    ij     ij 

If  the  interdiction  has  a  budget  of  K,  then  tlie  interdictor' s 

problem  is  to  find  a  set  of  m. .  which  minimizes  the  maximum 

flow  value  of  the  network  subject  to: 

E    C.  .  [u.  .  -  m.  .]  S  K 
(i,j)eA  iJ   iJ    iJ 

1.  .  ^m.  .  ^u.  .,  for  all  (i.j)  e  A, 
ij    ij    ij 

The  authors, like  Wollmer,  rely  on  the  original  network  being 
planar  and  undirected.   Also,  the  dual  network  is  used  in  this  method 
as  it  was- in  Wollmer' s.   A  very  brief  description  of  the  method  will 
now  be  given.   TJie  actual  statement  of  the  algorithm  is  rather  lengthy 
and  will  not  be  included. 
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The  first  step  of  the  algorithm  is  to  set  the  capacities  of 
each  arc  to  its  lower  bound  and  to  find  a  minimum  cut  for  this  "reduced" 
network  by  solving  a  shortest  path  problem  in  the  dual  network.   If 
this  solution  is  feasible  to  the  above  interdictor' s  problem,  then  it 
is  optimal ,  since  Diis  is  as  good  as  even  an  infinitely  wealthy  attacker 
could  do.   If  the  solution  is  not  feasible,  the  best  feasible  solution 
allowed  by  interdicting  only  the  arcs  of  the  minimum  cuL  is  found  and 
another  minimum  cut  to  the  reduced  problem  is  chosen  and  tested.   This 
procedure  ends  in  one  of  two  states:  either  a  minimum  cut  solution  is 
feasible  and  thus  is  optimal  or  no  minimum  cut  solution  is  feasible 
but  the  current  "best  solution"  has  been  noted.   If  the  latter  condi- 
tion holds,  the  cuts  of  the  reduced  network  are  tested  in  increasing 

th 
value  (actually  the  r    sliortest  path  problem  is  found  in  the  dual 

graph).   This  testing  and  updating  of  the  current  "best  solution"  is 

,,  th  ,, 

continued  until,  for  the  k    smallest  cut,   t)ie  best  possible  solution 

for  this  cut,  i.e.  ,  the  solution  witli  all  tlie  capacities  equal  to  the 

lower  bound,  is  not  better  than  the  current  best  feasible  solution. 

At  this  point,  the  algorithm  terminates  as  all  the  remaining  cuts  to 

the  reduced  problem  will  have  values  greater  than  a  feasible  solution 

of  another  cut. 

The  remainder  of  this  paper  is  divided  as  follows.   Chapter  2 

treats  the  single  most  vital  links  problem;  a  solution  method  for  this 

special  case  is  given  and  an  example  is  used  to  show  how  tlie  solution 

may  be  found.   Chapters  3  and  4  discuss  the  more  general  n  most  vital 
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links  problem.   Algorithms  are  developed  and  the  computational  experi- 
ence using  them  is  described  in  Chapter  5.   In  Appendix  A,  a  linear 
zero-one  and  a  zero-one  quadratic  formulation  of  the  n  most  vital  links 
problem  are  presented  and  discussed. 


CHAPTER  2 
THE  SINGLE  MOST  VITAL  LINK  PROBLEM 

Introduction 

The  problem  addressed  by  this  chapter  is  concerned  with  finding 
the  most  vital  link  in  a  single  commodity  flow  network  G(N,A)  (directed, 
undirected,  or  mixed).   An  arc  (x,y)e  A  is  declared  to  be  the  most 
vital  link  if  its  value  v(x,y)  is  at  least  as  large  as  the  value  of 
every  ot]ier  arc  in  the  network.   The  value  of  arc  (x,y)  is  defined  as 
the  difference  in  the  maximal  flow  values  in  networks  G(N,A)  and 
G(N,A  -  (x,y))  between  some  source  node  and  some  sink  node.   Thus,  v(x,jO 
reflects  the  reduction  in  the  maximal  flow  attainable  if  arc  (x,y)  is 
removed  from  the  network. 

Wollmer  has  developed  an  algorithm  [29]  for  determining  the 
most  vital  network  link.   This  algorithm  has  been  employed  by  Durbin 
[9]  to  determine  the  single  most  critical  link  in  a  highway  system. 

The  following  sections  of  this  paper  briefly  present  Wollmer' s 
method  for  finding  the  most  vital  link  in  a  network  and  develop  an 
improved  algorithm  for  solving  the  most  vital  link  problem.   An  example 
using  tlie  improved  method  is  included. 
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Wollmer's  Algorithm 

Wollmer's  algoritlim  for  finding  the  most  vital  link  in  a  net- 
work follows  as  a  consequence  of  the  following  theorem.  Tlie  proof  of 
this  theorem  is  given  in  reference  [29]. 

Tlieorem  (Wollmer): 

Suppose  the  network  G(N,A)  has  a  maximal  flow  value  of  v(F  ) 

while  the  network  G(N,A-  (x.y))  has  a  maximal  flow  value  of  v(F*  ). 

xy 

Then,  every  maximal  flow  pattern  of  G(N,A)  has  at  least  v(F  )  -v(F   ) 

xy 

units  of  flow  through  the  arc  (x,y).   Moreover,  there  is  a  maximal 
flow  pattern  of  G(N,A)  which  has  exactly  v(F  )  -v(F   )  units  of  flow 
over  the  arc  (x,y). 

As  Wollmer  points  out,  "the  above  theorem  reduces  the  problem 
to  one  of  finding  that  link  whose  minimal  flow  among  all  maximal  flow 
patterns  is  greatest."   Wollmer's  iterative  procedure  for  finding  this 
link  is  as  follows: 

STEP  0: 

Find  a  maximal  flow  pattern,  F*,  in  the  network  G(N,A)  and 
let  f*(x,y)  be  the  corresponding  flow  in  each  arc  (x,y)  e A, 

Set  the  "least  flow"  of  each  arc  (x,y)      equal  to  f  (x,y). 

Let  f*(p,q)  =  max  f*(x,y)  and  go  to  STEP  1. 
(x,y)  e  A 

STEP  1: 

Solve  a  maximal  flow  problem  for  the  network  G(N,A  -  (p,q)). 

Let  the  corresponding  maximal  flow  pattern  be  denoted  as  F 

'      "^  pq 

and  go  to  STEP  2. 
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STEP  2: 


(a)  Set  the  capacity  c(p,q)  of  arc  (p,q)  equal  to  v(F*)  -v(F*  ) 

and  solve  a  maximal  flow  problem  for  this  network  (i.e. ,  the 

network  G(N,A)  with  c(p,q)  =  v(F*)  -v(F*  )).   Call  the 

pq 

corresponding  maximal  flow  pattern  F  .   (Note  F   is  also 
a  maximal  flow  pattern  of  G(X,A).) 

(b)  Next,  compare  the  least  flov.-  of  each  arc  (x,y)  with  f  (x,y) 
and  if  f  (x,y)  <  least  flow  of  arc  (x,y) ,  replace  the  least 
flow  of  (x,y)  with  f  (x,y).   Reset  c(p,q)  to  its  original 
value  and  go  to  STEP  3. 


STEP  3: 


Let  U  =  max  [least  flow  of  arc  (x,y)}. 
(x,y)  eA 

(a)  If  U  S  v(F*)  -v(F*  )  terminate;  (p,q)  is  a  most  vital  link; 

(b)  If  U  >  v(F*)  -v(F*  );  find  an  arc  (p,q)  such  that  the 

pq 

least  flow  of  (p,q)  equals  U,  and  go  to  STEP  1. 


An  Improved  Algorithm 

Wollmer's  algorithm  considers  each  arc  as  a  candidate  for  the 
most  vital  link.   However,  a  necessary  condition  is  employed  in  the 
improved  algorithm  which  reduces  the  number  of  arcs  that  must  be  con- 
sidered ex-plicitly  as  candidates. 

Theoi'om  1: 

A  necessary  condition  for  an  arc  (a,b)  to  be  a  most  vital  link 
is  tliat  for  any  maximal  flow  pattern  in  t}ie  network  G(N,A)  ,  the  flow 
in  arc  (a,b)  is  at  least  as  great  as  the  flow  ovei'  every  arc  in  a  min- 
imal cut. 

Tlie  following  lenmia  will  be  useful  in  the  pi-oof  of  Theorem  1. 
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Lemma  2.1: 

If  (X,X)  is  a  minimal  cut  containing  at  least  two  arcs  in  a 
network  G(N.A)  and  if  arc  (x,y)  is  in  (X,X) ,  then  (X,X)  -  (x,y)  is 
a  minimal  cut  in  the  network  G(N,A-  (x,y)). 

Proof:   (Lemma  2.1) 

Suppose  that  (Y,Y)  is  a  minimal  cut  in  G(N,A-(x,y))  and  that 
C(Y,Y)  <C(X,X)  -c(x,y).   Note  that  (Y,Y)  U  (x,y)  has  to  be  a  discon- 
necting set  for  G(N,A)  and  that  C(Y,Y)  +  c(x,y)  <  C{X,X):   but  (X,X) 
is  a  minimal  cut  of  G(N,A)  and  (Y,Y)  U  (x,y)  is  a  disconnecting  set  of 
G(N,A).   Thus,  it  must  be  true  that  C(Y ,Y)  =  C(X,X)  -  c(x,y)  and,  hence, 
(X,X)  -  (x,y)  is  a  minimal  disconnecting  set  and  thus  a  minimal  cut  in 
G(N,A-  (x.y)). 

Proof:   (Theorem  2.1) 

Let  (X,X)  be  a  minimal  cut  in  G(N,A)  and  note  that  by  hypothesis 

arc  (a,b)  is  a  most  vital  link.   Assume  that  f*(a,b)<f  (p,q)=max  f  (x,y) 

(x,y)e(X,X) 
for  some  maximal  flow  pattern  F   defined  in  G(N,A).   It  will  be  shown 

that  this  assumption  leads  to  a  contradiction. 

(1)  Bj'  Lemma  2.1:   v(F*  )  =  v(F*) -f  *(p  ,  q)  ;  by  assumption  f  *(a  ,b)<f  *(p, 
therefore,  v(F*  )  =  v(F*)  -  f*(p,q)  <  v(F*)  -  f*(a,b). 

pq 

(2)  Further,  there  exists  a  flow  pattern  in  G(N,A  -  (a,b))  with  the 

value  v(F  ,)  =  v(F  )  -  f*(a,b),  and,  hence,  the  maximal  flow  value 
ab 

v(F*  ),  in  G(N,A-(a,b)) ,  must  satisfy:  v(F*  )  s v(F*)  -  f *(a,b) . 
ab  ab 

(3)  Conditions  (1)  and  (2)  implv  that  v(F*,  )>v(F*  );  but  this 

-  ab      pq 

leads  to  a  contradiction  since  by  assumption  (a,b)  is  a  most 
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vital  link  and  by  the  definition  of  a  most  vital  link  it 


follov.s  that:   v(F*  )  ^  v(F*  ),  for  all  (x,y)  e  A. 
ab       xy 


Q.  E.  D. 


Theorem  1  guarantees  that  those  arcs  whose  flow  is  less  than 
the  largest  flow  through  an  arc  in  a  minimal  cut  for  some  maximal  flow 
pattern  in  G(N,A)  need  not  be  considered  as  candidates  for  the  most 
vital  link. 

A  Labeling  Scheme 

A  labeling  scheme  is  employed  in  the  algorithm  for  finding 
a  most  vital  link  from  the  set  of  candidate  arcs.   At  the  outset  of 
the  labeling,  it  is  assumed  that  there  is  a  maximal  flow  pattern,  and 
minimal  cut  defined  for  G(N,A) .   Suppose  that  arc  (a,b)  is  an  arc  from 
the  candidate  set  (this  set  is  made  up  of  those  arcs  that  satisfy  the 
necessary  condition  of  Tlieorem  1).   Next,  an  attempt  is  made  to  label 
from  node  a  to  node  b  without  using  the  arc  (a,b).   Initially  node  a 
is  labeled  and  all  other  nodes  are  unlabeled.   The  labeling  scheme 
system.atically  searches  for  a  flow  path  from  node  a  to  node  b  which 
does  not  incl\ide  arc  (a,b)  such  that  flow  may  be  diverted  from  arc  (a,b) 
over  this  path.   If  node  b  is  labeled  (breakthrough),  such  a  path  exists. 
The  labeling  (and  backtracking)  rules  to  be  used  are  similar  to  those 
employed  by  Ford  and  Fulkerson  [11]  in  their  algorithm  for  the  solution 
of  the  maximal  flow  problem.   The  labeling  and  flow  changing  rules  are: 

Let  node  a  be  labeled  with  (-,  e  (a)  :^  ^') .      At  a  general  step, 
suppose  tlio  node  x  is  labeled  (z  ±  ,  e  (x))  and  tliat  the  nodes  x  and  y 
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are  connected  by  some  arc.   Tlien,  node  y  may  be  labeled  if  either  of  the 
following  situations  occur: 

(a)  f(y,x)  >  0;  then  node  y  is  labeled  (x  ,  e  (y))  where 
e  (y)  =  min(e  (x)  ,  f  (y,x))  ;  or 

(b)  f(x,y)  <  c(x,y);  then  node  y  is  labeled  (x  ,  e  (y))  where 
e  (y)  =  min  (  e  (x)  ,  c(x,y)  -  f  (x,y) ) . 

The  labeling  process  is  continued  until  either  node  b  is  labeled 
(breakthrough) ,  or  node  b  is  unlabeled  in  which  case  no  more  labeling  is 
possible  (non-breakthrough).   If  breakthrough  occurs,  node  b  must  have 

a  label  of  the  form  (q  ,  e  (b) )  for  some  node  q.   Likewise,  node  q  has 

±  -1- 

a  label  (r  ,  e  (q))  and  node  r  has  a  label  (p  ,  e  (r)),  etc.   Thus  a  series 

of  nodes  (and  a  path  of  arcs)  starting  with  node  b  and  ending  with  node 

a  is  defined. 

Let   e  =  min  [e  (b)  ,  f(a,b)].   For  each  arc  (x,y)  on  this  path, 

change  the  flow  as  follows: 

(a)  If  node  y  has  a  label  (x  ,  e  (y))  ,  replace  f (x,y)  by  f (x,y)  +  e; 

(b)  If  node  y  has  a  label  (x  ,  e  (y)),  replace  f(y,x)  by  f(y,x)  -  e. 
Finally,  decrease  the  flow  in  arc  (a,b)  by  e  units. 

As  Ford  and  Fulkerson  liave  shown,  the  labeling  schem.e  must  end 
in  one  of  two  mentioned  states:   breakthrough  is  achieved  or  breakthrough 
is  not  achieved.   The  following  results  indicate  what  can  be  deduced 
when  either  of  these  states  occur  at  the  end  of  the  labeling  process. 

Theorem  2: 

Let  F  be  a  maximal  flow  pattern  in  the  network  G(N,A).   Let  (a,b) 
be  an  arc  in  G(N,A)  and  use  the  labeling  scheme  to  label  from  node  a  to 
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node  b  without  using  arc  (a,b).   If  breakthrough  occurs,  the  value  of 
arc  (a,b)  is  less  than  or  equal  to  f  (a,b)  -  e.   If  non-breakthrough 
occurs,  the  value  of  arc  (a,b)  is  equal  to  f*(a,b). 

Proof: 

If  breakthrough  occurs  an  alternate  maximal  flow  pattern  is 

found  by  making  an  e  cliange  of  flow  in  the  path  established  by  the 

labeling  and  decreasing  the  flow  in  (a,b)  by  e.   Then, 

v(F*  )  >  v(F*)  -  f*(a,b)  +  e,  and,  since 
ab 

v(a,b)  ---    v(F*)  -  v(F*  )  ^  v(F*)  -  [v(F*)  -  f*(a,b)  +  e]  and 
ab 

e  >  0,  it  follows  that  v(a,b)  ^  f*(a,b)  -  e  <  f*(a,b). 

Assume  that  non-breakthrough  occurs,  then  by  the  nature  of  the 
labeling  rules,  no  maximal  flow  pattern  exists  in  G(N,A)  which  has  less 
than  f*(a,b)  units  of  flow  over  arc  (a,b)  and,  hence,  by  Wollmer's 
theorem  v(a,b)  =  f *(a ,b) , 

Q.  E.  D. 

Utilizing  these  results,  the  following  algorithm  can  be  employed 
to  locate  a  most  vital  link  in  a  network: 


STEP  0: 

(a)   Find  a  maximal  flow  pattern,  F  ,  in  tlie  network  G(N,A)  and  let 
(X,X)  be  a  minimal  cut.    Let  U*  =  max  c(x,y)  or,  alternately, 


(x,y)e(X,X) 
(x,y)e(X,X) 


*         * 
U   =  max  f  (x,y) 


The  Ford  and  Fulkerson  maximal  flow  algoritlim  [11]  may  be  used 
since  it  identifies  a  minimal  cut  as  well  as  a  maximal  Jlow  pattern. 
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(b)   Note  those  arcs  (x,y)  e  A  for  which  f  (x,y)  2:  U   and  store 
these  arcs  in  a  list;  these  arcs  form  the  candidate  set,  S. 
For  each  arc  in  this  set  define  an  upper  bound  as 
U(x,y)  :.  f*(x,y). 


STEP  1: 


Let  U(a,b)  =  max  U(x,y)  and  set  f(x,y)  -   f*(x,y)  for  all  arcs 

(x,y)  es 
in  G(N,A). 


STEP  2: 

Use  the  labeling  rules  to  label  from  node  a  to  node  b  without 
using  arc  (a,b). 

STEP  3: 

(a)  If  breakthrough  occurs,  use  the  backtracking  and  flow  changing 
rules,  replace  f(x,y)  by  the  resultant  flow  in  each  (x,y)  eA, 
and  if  f (a,b)  >  0,  repeat  STEP  2. 

(b)  Otherwise,  replace  U(a,b)  with  f (a,b)  and  if  U(a,b)  ^  max  U(x,y) 

(x,y)  es 

terminate;  arc  (a,b)is  a  most  vital  link. 

Otherwise,  replace  f  (x,y)  with  f (x,y)  for  all  (x,y)  e  A 

and  go  to  STEP  1. 

By  the  use  of  this  algoi'ithm,  a  maximal  flow  pattern  is  always 
maintained  and  once  non-breakthrough  results  for  the  arc  (a,b)  and 


An  alternative  to  STEP  3a  is  to  test  if  f (x,y) <U(x,y) , 
(x,y)  es,  and  if  f(x,y)  ^U*,  to  replace  U(x,y)  with  f(x,y);  if 
f (x,y)  <U*,  then  arc  (x,y)  may  be  dropped  from  the  set  S.   The  com- 
putations are  continued  by  either  repeating  STEP  2  if  arc  (a,b)  has 
not  been  dropped  or  returning  to  STEP  1  if  arc  (a,b)  has  been  dropped. 
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U(a,b)  Si  U(x,y)  for  all  arcs  (x,y)  of  the  candidate  set  S,  the  arc 
(a,b)  can  be  declared  a  most  vital  link,  i.e.  ,  non-breaktlirough  implies 
that  v(a,b)  =  f*(a,b)  =  U(a,b)  and,  thus,  v (a ,b)  >U(x,y)  > v(x,y)  for 
(x,y)  e  S  and  by  Theorem  1  tlie  most  vital  link  of  G(N,A)  is  an  element 
of  S. 

The  other  aspects  of  the  method  tliat  must  be  considered  are: 
(1)  is  the  method  finite  and  (2)  is  there  a  way  to  locate  alternative 
optimal  solutions?   The  finiteness  of  the  procedure  follows  since  there 
can  be  only  a  finite  number  of  arcs  in  the  candidate  set  and  there  is 
a  finite  number  of  label ings  that  can  be  made  for  each  arc  of  the  set. 
Alternative  optima   are  readily  identified  as  follows:   after  an  optimal 
solution  has  been  found  and  other  candidate  arcs  exist,  reapply  the 
algorithm  to  any  arc  in  tlie  candidate  set  whose  upper  bound  is  equal  to 
the  value  of  the  most  vital  link  (a,b),  i.e. ,  set  U  =  U(a,b),  delete 
arc  (a,b)  from  the  candidate  set,  and  reapply  tlae  algorithm.   If  upon 
reapplication  of  the  algoritlim  the  most  vital  link  has  a  value  equal 
to  U  ,  then  an  alternative  solution  has  been  found.   In  the  latter  case, 
the  algorithm  is  reapplied  using  a  further  reduced  candidate  set  and 
the  procedure  is  continued  until  all  alternative  optimal  solutions 
have  been  four.d. 


An  Example 


Consider  the  flow  network  G(N,A)  shown  in  Figure  2. 1. 


23 


10,10 


Figure  2.1,   Maximal  Flow  in  Sample  Network  2.1  with 
Source  Node  s  and  Sink  Node  t. 


Associated  with  each  arc  in  the  network  is  an  ordered  pair  of 
numbers.   The  first  number  corresponds  to  the  capacity  of  the  arc  and 
the  second  number  corresponds  to  the  amount  of  flow  over  the  arc  in 
some  maximal  flow  pattern. 

The  improved  algorithm  will  be  employed  to  find  all  most  vital 
links  in  the  network  G(N,A). 
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STEP  0: 

(a)  Figure  2. 1  presents  a  maximal  flow  pattern  F  in  the  network 
with  27  units  of  flow  from  node  s  to  node  t.  The  individual 
arc  flow^s  are: 

f*(s,a)  =  25  f*(a,c)  =  15 

f*(s, b)  =   2  f*(a,d)  =  10 

f*(c,t)  =  17  f*(b,c)  =   2 

f*(d,t)  =  10  f*(b,d)  =   0 

The  corresponding  minimal  cut  (which  is  unique)  contains  the 
arcs  (a,c),  (b,c),  and  (d,t).   Hence, 

U*  =  max  f*(x,y)  =  f*(a,c)  =  15 
(x,y)e(X,X) 

(b)  The  set  of  candidate  arcs  is: 

S  =  {(s,a) ,  (c,t) ,  (a,c)} 

STEP  1: 

Start  with  arc  (s,a),  since  U(s,a)  -   f  (s,a)  -   max  U(x,y). 

(x,y)  es 

STEP  2: 

Labeling  from  node  s  to  node  a  without  using  arc  (s,a)  is 
possible  over  the  path  containing  the  arcs  (s,b),  (b,d), 
and  (a,d).   The  value  ot    e  -    8. 

STEP  3: 

(a)   Since  breaktlirough  occuri-ed,  the  e  units  of  flow  are  removed 
from  arc  (s,a)  and  added  to  each  arc  in  tlie  path  found  in 
S'J'EP  2.   The  revised  maximal  flow  jjattern  is  shown  in  Figure  2.2. 


25 


STEP  2; 


Figure  2.2.   Revised  Maximal  Flow  in  Sample  Network  2.1, 


Starting  with  the  maximal  flow  pattern  of  Figure  2.2,  it  is 
not  possible  to  label  from  node  s  to  node  a  and  non- 
breaktlirough  has  occurred. 


STEP  3: 


(b)   Replace  U(s,a)  with  the  flow  f(s,a)  =  17.   Since 

U(s,a)  '^.   max  U(x,y)  =  U(c,t),  arc  (s,a)  is  a 
(x,y)  e  S 

most  vital  link. 
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Since  U(c,t)  =  U(s,a)  =  17,  arc  (c,t)  may  also  be  a  most  vital 
link.   In  order  to  test  this  hypothesis,  arc(s,a)  is  dropped  from  the 
set  S  and  the  algoritlim  is  reapplied  starting  witli  the  maximal  flow 
pattern  in  Figure  2.2. 

Employment  of  the  labeling  rules  from  node  c  to  node  t  results 
in  non-breakthrough  and  arc  (c,t)  is  an  alternative  most  vital  link. 


CHAl'TER  3 
THE  n  MOST  VITAL  LINKS  PROBLEM 

Introduction 

As  stated  in  CTiapter  1,  determining  the  n  most  vital  links  of 
a  network  involves  locating  that  set  of  n  arcs  whose  removal  most 
impairs  the  throughput  capability  of  the  remaining  network.   In 
Chapter  2,  results  pertaining  to  the  special  case  where  n =  1  were 
given.   While  the  algorithm  developed  in  that  chapter  is,  with  minor 
changes,  valid  for  the  more  general  problem  it  is  not  computationally 
feasible  for  large  problems  with   n>l.   Furthermore,  as  Wollmer  in  [29] 
notes,  the  single  most  vital  link  may  not  be  one  of  the  n  most  vital 
links.   Thus,  the  single  most  vital  link  problem  will  be  treated  as 
a  special  case  and  another  approach  will  be  employed  for  the  more  gen- 
eral problem. 

A  Cut  Problem  Approach 

In  this  chapter  the  n  most  vital  links  problem  is  formulated  as 
a  cut  problem  and  an  algorithm,  which  uses  the  maximum  flow  algorithm 
of  Ford  and  Fulkerson  [11], is  developed  to  solve  this  problem. 

Theorem  3.1  and  Corollary  3,1.1  characterize  the  n  most  vital 
links  of  G(N,A)  as  elements  in  a  cut  between  the  source  and  sink  of  G(N,A), 
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Theorem  3.1: 


If  the  arcs  in  A   are  the  n  most  vital  links  in  the  network, 

G(N,A) ,  then  there  exists  a  cut  a'  such  that  a'  e  a'  and  a'-a'  is 

n  -  n 

a  minimum  cut  in  G(N,A-A  ). 


Proof: 

Find  the  minimum  cut  in  G(N,A-a') ,  say  a"  with  capacity  CCa") 

n        n  n 

Now  consider  the  arc  set  A."   U  A  ,   This  is  certainly  a  disconnecting; 

n    n  ° 

set  for  G(N,A),   Suppose  it  is  not  a  cut.   Then  there  exists  at  least 

one  arc  a   e  a''  U  a'  such  that  k"   \J  k'    -   fa.]  is  also  a  disconnecting 
inn  nn]-"  "= 

set  for  G(N,A).   If  a.  e  A  ,  then  a"  is  not  a  minimum  cut  for 
in        n 

G(N,A-a').   If  a   e  a',  then  a'  U  [a,  ]  -  fa.],  where  a,  is  any  arc 
n        in        n     ki"^         K 
// 

m  A^,  IS  a  better  set  of  n  links  to  remove  from  G(N,A).   In  either 


case  a  contradiction  results.   Hence,  by  letting  a'  =  a'  U  k" ,    the 

n    n 

result  follows. 

Q.  E.  D. 
Corollary  3.1.1: 

C(A  )  >  C(A*)  for  all  arc  sets,  A*  c  a' ,  which  contain  at  most 
n       n  n  ~ 

n  arcs. 

Proof: 

Suppose   C(A^'^)  >  C(a')  for  some  A*  c  a'  containing  at  most   n 
arcs.   Note  that  for  any  A*  c  a' 

C(a')  =  C(A*)  +  C(A'-A*)  =  C(a')  +  c(a'-a'), 
n  n       n         n 

hence,  CCa'-a"^)  <  C(a'-a'),  but  this  implies  that  a'  is  not  the  n  most 
n  n  n 

vital  Ijnks  in  G(N,A)  which  is  a  contradiction. 

Q.  E.  D, 
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Corollary  3. 1.2: 

C(A  )  ^   C(A   )  where  A    is  any  minimum  cut  of  G(N,A)  and  a' 
n       n  r 

/// 
contains  the  n  largest  arcs  of  A   .  . 


Proof: 


This  result  follows  since  C(a'-a')  ^   C(a'"-a'") 


(or  C(A  )  -  C(A  )  ^  C(A   )  -  C(A   ))  and   C(A  )  S  C(A   ). 
n  n 

Q.  E.  D. 

Note  that  due  to  Theorem  3.1  and  Corollary  3.1.1,  the  n  most 

vital  links  problem  of  G(N,A)  can  be  stated  as  the  following  cut  problem: 

Find  the  n  largest  arcs,  A*,  of  the  cut  A*  which  satisfies: 

n 

C(A*-A*)  =  min  C(a'-a') 
n  n 

for  all  cuts  A 

where  A   denotes  the  n  largest  arcs  of  A  . 
n 

Development  of  a  Cut  Algorithm 
Let  G  (N,A)  be  the  network  G(N,A)  with  the  capacity  of  each 


arc  a.  e  A  defined  as:        , 
1  / 


c  (a.)  = 


c(a . )     if   c(a. )  <  u 


if   c(a  . )  S:  u. 


Note  that  G  (N,A)  s  G(N,A) .   For  any   \  ^A   let 


C  (A  )  =  l^      c"(a.). 

a.  eA, 
1   k 
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Theorem  3.2: 

If  A   is  a  minimum  cut  in  G  (N,A)  with 
a'  s  fa   a  eA   and   c  (a  )  =  u],   and   a  =   A    is  defined  as  the 

U       1     1  1  u 

number  of  arcs  contained  in  A  ,  then  the  arcs  in  A   are  the  a   most 

u  u 

CO 

vital  links  in  G  (N,A). 

Proof: 

A   contains  the  a   largest  arcs  in  A   for  both  G  (N,A)  and 

CO  y  CO 

G  (N,A).   Suppose  A   does  not  contain  the  a   most  vital  arcs  in  G  (N,A) 
Then,  from  Theorem  3.1  there  exists  a  cut  A*  such  that: 

CO      -(;       i  CO  f  , 

C  (A  -A  )  <  C  (A  -A  )  (1) 

where  A   are  the  a   largest  arcs  of  A*  in  G  (N.A). 
u 

CO  f  / 

Note  that   c  (a.)  <  u   for  all   a.  eA  -A  ,  since 


c  (a.)  >  u  ^  c  (a.)  =  u  =>  a  eA  « 

1  1  1    u 

Hence, 

^   /   /      u   /   / 
C  (A  -A  )  ^  C  (A  -A  )  .  (2) 


Equations  (1)  and  (2)  imply  that 


but 


Hence, 


"*   *   *      u   /   / 
C  (A  -A  )  <  C  (A  -A  ) 


U      *      +  CD     s|c      ;); 

C  (A  -A  )  S  c  (A  -A  ) 


u   *   *      u   /   / 
C  (A  -A  )  <  C  (A  -A  ) 


C"(A*)  -  C^(A*)  <  c"(a')  -  c"(a') 
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Now,  since  c  (a.)  ^  u  for  all  a  G A*  and  c  (a  )  -   u, 
1  1    u         i 

for  all   a.  e  a'  ,    c"(A*)  ^  c"(a')  whicli  implieE  that   c"(A*)  <  c"(a') 
1    u      u        u 

which  is  a  contradiction,  since  A   is  a  minimum  cut  of  G  (N,A). 

Q.  E.  D. 

Corollary  3.2. 1: 

The  maximum  flow  through  the  network  G(N,A-a')  is  equal  to 

cca'-a'). 

u 

Proof: 

A   -  A   is  a  minimum  cut  in  G  (N,A-A')  and,  since 
u  u 

C  (a'-A')  =  C  (a'-a')  ,  k'  -  a'      is  also  a  minimum  cut  in  G  (N,A-A'). 
u  u         u  u 

Q.  E.  D. 

Corollary  3.2.2: 

Let   a'  =  fa.   a.  sa'   and   c  (a.)  >  u]   where   A '   =  a 
X      1    1  1      ^  '  x' 

and  a'  =  [a.   a.  ga'   and   c  (a.)  =  u]   where   a'   =  3  .   Then 
p    ^  1    1  1      •"  '  p' 

k 
A   U  a.    are  the  Ck*  +  k  most  vital  links  of  G(N,A)  for   a  eA  and 
'^  i.l  ^  '         X    p 

any  k  =  0,1,2,. ..,3. 

Proof : 

Same  as  for  Theorem  3.2. 

Q.  E.  D. 

Let  a  minimum  cut  of  G  (N,A)  be  denoted  by  A   and  define  A 

'  ^  u 

to   be   the   set   of    arcs   defined  by     A      ={a|aeA        and      c    (a)    =   u]. 
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Theorem  3.3: 

Let  T\   be  defined  as  the  smallest  number  of  arcs  contained  in 
any  cut  (i.e.  ,  T]  is  the  minimum  number  of  arcs  necessary  to  disconnect 
the  network  G  (N,A)),  then 

T]  -    \\l    ^  l\l    ^  ...    s   \aJ    ^  ...    ^   |aJ    ^  0. 

Proof: 

From  the  definition  of  A   it  follows  that   A   s  0.   Also, 

U  'col  ' 

by  noting  that  in  G  (N,A)  all  arcs  have  a  capacity  of  1,  it  follows 
that  the  minimum  cut  of  G  (N,A)  must  be  tliat  cut  which  contains  the 


fewest  arcs  (i.e.,   T]  =  [a  |). 


It  must  now  be  shown  that   p<q  implies  |A  |  ~i:    (a  ).   Assume 

that  this  is  not  the  case,  i.e.  ,  assume  that   p<q  and   A   <  A   . 

p'    '  q 

Note  that  by  Theorems  3.1  and  3.2  the  arcs  in  A   are  the  A 

P  P' 

most  vital  links  of  G(N,A)  and  the  maximal  flow  value  of  G(N,A--A  ) 

P 

is  C(A  -A  ).   Similarly,  the  arcs  in  A   are  the  A   most  vital  links 
P  q  q' 

of  G(N,A)  and  removing  A   from  G(N,A)  leaves  a  maximal  flow  value  of 

q 

C(A^-A  )  in  G(N,a"'-A  ).   Also,  since  it  is  assumed  that  Ia  I  >  Ia  I  , 
q  q  '  q'    '  p'  ' 

it  follows  that 

CO    r)  OO    Q 

C  (A  -A  )  :>  C  (A^-A  )  ,  (3) 

p         q 

Since  for  any  arc   a  eA  ,  c  (a)  >-   p  if  and  only  if   a  eA  , 


it  follows  that 


c"'(a''-A  )  =  d^(A'^-A  )  , 
P  P 


(where   d'(A''-A  )  =  d'(A^)  -  p|a  |), 
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Similarly,  it  follows  that 


C°°(A^-A  )  --  C^(A^-A  )  =  C*^(A^)  -  cJa  | 
q  q  '  q' 


Using  these  identities  the  inequality  (3)  can  be  written  as 


cP(aP)  -  p|a  U  C^A^)  -  q|A  |  .  (4) 

p  q 


Finally,  note  that 


or  equivalently 


C^(A^-A  )  S:  C^(A^-A  )  , 

q        q 


C^A^)  -  q|Aj  S  CP(A^)  -  pIaJ  . 

Thus,  using  this  last  relationship  with  (4)  and  the  assumption  that 

\a    I  >  |a  I ,  it  follows  that 
'  q'     P 

cP(aP)  >  cP(A^)  . 
But  this  is  a  contradiction,  since  A   is  a  minimum  cut  of  G  (N,A) 


and,  thus,  Ia  I  ^  |a  | 


Q.  E.  D. 


The  Modified  Cut  Method 

Theorems  3,2  and  3.3  provide  the  basis  for  a  very  simple 
algorithm  to  solve  the  n  most  vital  links  problem  for  some  values 
of  n.   The  algorithm  consists  of  defining  a  graph  as  in  Theorem  3.2 
and  then  systematically  changing  the  value  of  u.   Stated  formally, 
the  algorithm,  "the  modified  cut  method,"  becomes: 
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STEP  0: 

Let   II  =  1   and  go  to  STEP  1. 

STEP  1: 


■  Set   c  (a.)  = 


c  (a. )     if   c  (a. )  <  u 


if   c(a.)^u 


Solv'e  for  the  minimum  cut  of  G  (N,A)  . 
Note  a",  a  ,  and  Ia  I .   Go  to  STEP  2. 

STEP  2: 

If   |a  I  >  n   let   u=  u  +  1  .   Go  to  STEP  1. 
u 

If   A  I  =  n   terminate  as  tlie  n  most  vital  links  are  the 
'  u' 

arcs  contained  in  A  . 
u 

If   A   <  n   terminate,  a  gap  has  occurred  at  the  n  most 
'  u' 

vital  links  problem. 


"^Note  that,  due  to  Corollary  3.2.2,  if  k  of  the  |a  |  ceiling 

arcs  have   c  (a  )  =  u,   then  the  A  I  -  k   through  the  A   most  vital 
i  u  u' 

links  have  been  found.   Tliis  addition  to  the  algorithm  has  proved 
useful  in  the  examples  that  have  been  tried  to  date. 
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If  the  algorithm  terminates  with  (a  I  ::;  n,  then  from  Theorem  3.2 

u' 

the  n  most  vital  links  have  been  found.   If  the  algorithm  terminates 

with   |A  I  <  n,  then  there  is  a  "gap"  at  the  value  of  interest  (i.e., 

the  algorithm  obtains  some  sets  of  most  vital  links  but  not  the  n  most 

vital).   This  problem  will  be  treated  in  detail  in  Chapter  4.   The 

algorithm  clearly  terminates  in  a  finite  number  of  steps,  since  for 

u  =  max  c(a.)  +  1  for  all   a.  eA,  |a  I  =  0  (for  every  minimum  cut  of 
1  1        u 

G  (N,A)).   It  should  be  noted  that  for  u  =  1  all  arcs  in  A   are  at 

capacity  u  =  1.   Hence,  at  least  for  this  value  the  algorithm  finds 

a  set  of  most  vital  links.   In  particular,  it  finds  the  fewest  arcs 

necessary  to  disconnect  the  network.   Also,  from  Theorem  3.3,  \a        is 

u 

non-increasing;  hence,  if  |A  |  <  n,  there  is  no  point  in  investigating 
a  larger  value  of  u. 

Finally,  observe  that  the  maximum  flow  pattern  found  for  some 
G  (N,A)  is  feasible  for  G    (N,A)  and  can  be  used  as  an  initial  flow 
in  finding  the  maximum  flow  in  G    (N,A). 

An  Example  Using  the  Modified  Cut  Method 

The  modified  cut  method  will  be  illustrated  by  determining  the 
two  most  vital  links  in  Examjjle  Network  3.1  as  defined  in  Figure  3.1. 
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Figure  3.1.   Example  Network  3.1. 


Note  that  G  (N,A)  is  defined  as  in  Figure  3.2. 


Figure  3.2.   G  (N,A)  of  Example  Network  3.1. 
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'Hie  minimum  cut  value  of  G  (N,A)  is  4  with   A  ^  f  a,  .  ^o  .  ^r, .  ^  ^  ] 

1   z   J   4  -^ 

a  minimum  cut  of  G  (N,A)   and   |A  |  =  4. 

2 
For   G  (N,A)  (see  Figure  3.3),  the  desired  minimum  cut  is 

A  =  {'i,  1  ' ''-,  o  ' ''i -? ' ''i /i}  ^^^      l^nl    =  ^-   The  minimum  cut  value  of 

11    xA         1*3    I'l  ^ 

G  (N.A)  is  7. 


Figure  3.3.   G  (N,A)  of  Example  Network  3.1, 


For   G  (N,A) ,  as  shown  in  Figure  3.4,  the  desired  minimum  cut 

is  again  A  =  {'^-i  i  >  ^io  '  ^i -5 '^i/i]  '   l"^-:!'  "  ^   ^"^  *^®  minimum  cut  value 
is  9.   Thus,  a    and  a    are  the  two  most  vital  links  and  the  problem 

is  solved.   It  may  be  noted  that  if  the  result  of  Corollary  3,2.2 

2 
were  invoked  at  G  (N,A) ,  the  same  conclusions  could  liave  been  made, 

since  a  ^  was  in  the  minimum  cut  of  G  (N,A)  ;  |A  j  =  3;  and,  c  (s.   ^)   =2. 
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Figure  3.4.   G  (N,A)  of   Example  Network  3.1. 


An  Example  V>"nere  a  Gap  Occurs 


In  stating  the  modified  cut  method  it  was  noted  that  the  algo- 
rithm may  end  in  a  gap,  that  is,  jA  |  >  a  and   |A  +1|  <  o^. 
In  Chapter  5  the  frequency  of  occurrence  of  these  gaps  in  the  sample 
problems  to  date  is  noted  and  discussed.   In  this  section  a  network 
will  be  used  for  which  the  modified  cut  method  fails  to  find  the  n  most 
vital  links. 

Consider  the  Example  Network  3.2  as  defined  in  Figure  3.5  and 
use  the  modified  cut  method  to  find  the   two  most  vital  links  of  that 
network. 
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Figure  3.5,   Example  Network  3.2. 

Note  that  for  G  (N,A)  the  arcs  leaving  node  s  form  the  minimum  cut. 
(See  Figure  3.6), 


Figure  3.6.   G  (N,A)  of  Example  Network  3.2. 
Thus,  tlie  minimum  cut  value  of  G  (N,A)  is  3  a}id  |A  J  -    3. 
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G  (N,A)  is  noted  in  Figure  3.7.   Again  the  arcs  leaving  node  s 
form  the  minimum  cut,  [A  |  -  3,  and  the  minimum  cut  value  is  6. 


Figure  3.7,   G  (N,A)  of  Example  Network  3.2. 


The  arcs  that  enter  node  t  form  a  minimum  cut  for  G  (N,A) 
as  shown  in  Figure  3.8.   This  minimum  cut  has  a  cut  value  of  8, 
Each  of  the  arcs  in  this  cut  has  a  capacity  of  2  in  the  original 
network,  so   |a  |  =0,  and  a  gap  exists  at  the  desired  problem. 
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Figure  3.8.   G  (X,A)  of  Example  Network  3.2. 


In  Chapter  4  the  gap  problem  wi] 1  be  treated  and  the  resulting 
algorithm  will  be  applied  to  Example  Network  3.2  given  in  this  chapter. 


CHAPTER   4 
RESOLUTION  OF  THE  GAP  PROBLEM 

Introduction 

In  this  chapter  the  gap  problem  introduced  in  Chapter  3  is 
addressed.   As  shown  in  Example  3.2,  at  termination  the  modified  cut 
method  may  not  determine  the  desired  n  most  vital  links,  i.e.  ,  the 
method  may  end  in  a  gap.   Results  will  now  be  developed  and  an  algo- 
rithm presented  wliich  will  resolve  a  gap  if  it  occurs. 

It  will  be  assumed  in  this  chapter  that  the  arcs  of  a  network 

are  ordered  such  that  if  a.  and  a.  are  arbitrary  arcs  of  G(N,A)  and 

if   i  <  j,  then   c(a.)  ^   c(a.).   This  leads  to  no  loss  of  generality, 

since,  if  the  condition  does  not  hold  in  the  original  network,  the 

arcs  can  be  renumbered  without  changing  the  network.   Also,  since  in 

what  follows  only  disconnecting  sets  between  the  source  node  s  and 

sink  node  t  are  of  interest,  the  term  disconnecting  set  will  imply  a 

disconnecting  set  between  s  and  t.   If   F  is  a  disconnecting  set, 

then  C(F-F  )  will  be  called  the  n-reduced  value  of  F,  where  F 
n  n 

denotes  the  n  largest  arcs  of  F.   Finally,  if   C(F-F  )SC(H-H  )  for 

n        n 

all  disconnecting  sets  in  some  set  J/,  then  F  will  be  called  the 
weakest  n-disconnecting  set  of  %C. 
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Some  Preliminary  Results 

Recall  that  in  Chapter  3  the  n  most  vital  links  were  shown  to 

be  the  n  largest  arcs  A*  of  a  disconnecting  set   A   which  satisfies: 

C(A  -A  )  ^  C(a'-a')  for  all  disconnecting  sets   a'  in  Ci(N,A).   Thus, 
n  n 

solving  for  the  optimal  solution  of  this  disconnecting  set  problem 
determines  the  n  most  vital  links.   Initially,  every  disconnecting  set 
of   G(N,A)  is  a  candidate  for  the  optimal  solution  to  the  most  vital 
links  cut  problem.   The  following  theorem  can  be  used  to  eliminate  from 
further  consideration  many  of  the  disconnecting  sets  of  G(N,A), 

Theorem  4.1: 

Let   a  ,  a   be  arcs  of  G(N,A)  such  that   c(a.)  >  c(a.),  and 
i    J  1       J 

define  the  set  U   to  be: 

J^  s  [h  I  H  is  a  disconnecting  set  of  G(N,A)  and  a.  e  H  ]. 

If   H   is  a  disconnecting  set  of  U   which  satisfies: 

C(H  -H  )  ^  C(H-H  )  for  all  disconnecting  sets  U    e  U 
n         n 

and  if  there  exists  a  disconnecting  set  F  such  that   a.  e  F   and 

C(F-F  )  <  C(H  -H  ),  then   a.  0   F, 
n  n  1 


Actually  only  cuts  were  considered  in  Chapter  3,  but  since 
the  weakest  n-disconnecting  set  of  G(N,A)  is  a  cut  and  every  cut  is 
a  disconnecting  set,  this  generalization  follows. 


44 


Proof: 

Given   that      a      e  F    ,    there   are   three  possible   cases 

Case   1   —      a.    €  F        and      a.    e  F 
J  n  in 

Case  2   —      a      e  F        and      a.    e  F-F 
J  n  in 

Case  3   —      a.    e   F        and      a,    ft  Y. 
J  n  1 

If      Case    1  holds,    then   F    e  ?^,    but   by    the  hypothesis      C(H   -H   )  ^  C(H-H   ) 

for   all      He?/     so   this    case   cannot   occur.       Case   2    is   not   possible, 

since,    if      c(a    )    >   c(a    ),    then      a      e  F-F        implies   that      a.    e  F-F    , 
1  J  in  J  n 

and   if      c(a)=c(a),    then     F    €  J/     and   the   comments   made   about   Case    1 
i  J 

are  valid   and   lead   to   a   contradiction  for   this    case.      Thus,    only 
Case   3    is   possible   and   the  result   follows. 

Q.    E.    D. 


Corollary   4.1.1: 

Given  the  network  G(N,A),  let  Y  =    fa  , a„ a  ]  be  a  set  of 

^12      p-* 

arcs  such  that  F  c:  A  and  let   a,  e  A  where   k>p  (i.e.,  c(a  )  ^  c(a  ) 
~  k  k       1 

for   i  =  l,2,.,.,p).   If  among  all  disconnecting  sets  H  which  satisfy 

H   n  r  *'  5 ,  H   is  the  weakest  n-disconnecting  set  and  if  there  exists 
n 

a  disconnecting  set  F  for  which  a,  e  F   and  C(F-F  )  <  C(H  -H  ) ,  then 

k    n  n  n 

F  n  r  =  $. 


Proof: 

Since   c(a.)  S:  c(a.  )   for   i  =  1,2 p,  the  result  follows 

1       k 

from  Theorem  4.  1. 
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Let  "U   and  ^  be  sets  defined  by: 

%C    =  {}I  I  H  is  a  disconnecting  set  of  G(N,A)  and  a   €  H  } 

*   r  *  I  *  *•> 

V  =  [11  |H   is  a  disconnecting  set  of  G(N,A)  and  a   «!  H  ]. 

If  the  weakest  n-disconnecting  sets  in  both  U  and  ?/  are  found,  then, 
by  Theorem  4.1,  the  n  most  vital  links  of  G(N,A)  are  determined.  The 
problem  of  determining  these  arc  sets  will  now  be  investigated. 

The  Extended  Modified  Cut  Method 

Observe  that  in  defining  the  modified  cut  method  it  was  only 
necessary  to  find  the  weakest  n-cut  of  G(N,A) ,  i.e. ,  only  cuts  of 
G(N,A)  and  not  all  disconnecting  sets  were  considered.   In  this  section 
an  extended  version  of  this  method  will  be  employed,  and  it  will  be 
necessary  to  consider  disconnecting  sets  rather  than  cuts  of  G(N,A). 
That  is,  for  a  fixed  set  of  arcs  A  the  n-|A(  most  vital  links  will  be 
found  in  G(N,A-A).   Thus,  while  the  weakest  (n-j a| ) -disconnecting  set 
of  G(N,A-A)  will  be  a  cut  of  G(N,A-A),  the  set  together  with  A  may  only 
be  a  disconnecting  set  of  G(N,A).   This  will  cause  no  difficulty,  how- 
ever, as  all  the  results  which  have  been  stated  arc  valid  for  both 
disconnecting  sets  and  cuts. 

Observe  that  the  weakest  n-disconnecting  set  of  the  set  ^  may 
be  characterized  as  the  set  whose  n  largest  arcs  consist  of  a   and  the 
n-l  best  arcs  to  remove  from  G(N,A) ,  given  a   is  to  be  removed.   But, 
this  is  exactly  the  disconnecting  set  which  would  be  found  by  setting 
A  =  a   and  solving  for  the  n-(A|  most  vital  links  of  G(N,A-A).   Thus, 
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the  modified  cut  method  can  be  used  to  find  the  weakest  n-disconnecting 
set  of  J/.   Assuming  that  the  desired  set  of  JJ^  is  found  (the  case  where 
a  gap  occurs  in  this  problem  will  be  treated  later  in  this  section) , 
the  set  J/  will  now  be  considered,  and  it  will  be  demonstrated  that  the 
modified  cut  metliod  can  also  be  extended  to  find  the  weakest 
n-disconnecting  set  among  those  sets  which  do  not  contain  arcs  from 
some  fixed  set  P.   (Thus,  for  F  =  a  ,  the  method  to  be  presented  will 
locate  the  weakest  n-disconnecting  set  of  )i/'  .  ) 

It  is  assumed  in  what  follows  that  there  exist  disconnecting 

sets  in  G(N,A)  which  do  not  contain  arcs  from  F,  since,  if  this  is  not 

*    * 
the  case,  then  in  the  context  of  the  last  section  /!/     =    9 . 

Given  the  fixed  set  of  arcs  F  and  G(N,A),  let  ^  (N,A)  be  G(N,A) 

u 
with  a  capacity  function  C   defined  by: 


CO    if   a.  e  F 

1 

C(a)-(    u    if   a   jSF   and   c(a.)  >  u 
i     1  1  1 


,c(a.)   if   a.  fi  r      and   c(a.)  <  u 
11  1 


OyU 

An  arc  a.  will  be  called  a  u-ceiling  arc  of  c(a.)  =  u  in  aJ  (N,A), 


Theorem  4.2: 


If  a  minimum  cut  A  of  £>     (N,A)  contains  oi   u-ceiling  arcs, 
then  C(A  -A  )  ^  C(B-B  )  for  all  cuts  B  such  that  B  fl  F  =  $. 
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Proof: 


*  n  U  U  *  11 

A      is    a  minimvim   cut   of    afr     (N,A),    therefore  C    (A   )    -^   C    (B)    for 


all    cuts    B.      Also,    for   all    cuts   B   that    satisfy      B  H  F  =    $,    C    (B    )    ^   ua. 
But,    by    the  hypothesis   C    (A    )    =   ua;    therefore, 

u**  u*  ^u  u 

C    (A    -A    )    =    C    (A    )    -   uoi   ^  C    (B)    -   ua  S    C    (B-B    )  (1) 

Q-  a 


for  all  cuts  B  such  that   B  H  F  =  $. 

Next,  observe  that  for  cuts  B  which  satisfy   B  R  F  =  $, 

u  * 

C(B-B  )  ^  C  (B-B  )   and  that,  since  A   contains  exactly  O"  u-ceilinpr 

a  a 

*   *     u   *   * 
arcs,  C(A  -A  )  =  C  (A  -A  ) .   Combining  these  last  conclusions  with  (1) 

gives  the  desired  result. 

Q.  E,  D. 


Recalling  Theorem  3.3,  which  is  valid  for  the  «^  (N,A)  as 
well  as  the   G  (N,A)  networks   and  Theorem  4.2,  the  extended  modified 
cut  method  may  be  stated  as: 
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STEP  0: 


Note  the  fixed  set  P. 

Let   u  =  1   and  go  to   STEP  1. 


STEP  1: 


Set   c  (a .  ) 


if   a.  e  r 

1 

CO  CO 

c  (a  )   if   c  (a.)  <  u,  and  a .    fL  V 
1  11 

CO 

u     if   c(a.)Su,  and   a.  ft  i 


STEP  2: 


Solve  for  the  minimum  cut  of  it    (N,A) ;  if  the  minimum  cut  value 

is  not  finite  stop  (all  cuts  contain  arcs  of  D  ,  otherwise 

note  a",  A  ,  and  Ia  1.   Go  to  STEP  2. 
u        u 

If   Ia  I  >  n   let   u  =  u+  1.  Go  to  STEP  1. 
'  u' 

If   A   =  n   terminate  as  the  n  most  vital  links  are  the 
u' 

arcs  contained  in  A  . 
u 

If   [a  I  <  n   terminate,  a  gap  has  occurred  at  the  n  most 
vital  links  problem. 


Finally,  the  problem  of  resolving  gaps  that  arise  in  using  the 
extended  modified  cut  method  to  find  the  weakest  n-disconnecting  set 
in  either  ^  or  %(     must  be  addressed.   Suppose  a  gap  occurs  in  solving 
for  the  n-jAj  most  vital  links  in  G(N,A-A)  where  the  disconnecting  sets 
containing  arcs  from  some  fixed  set  Y   are  not  to  be  considered.   This 
is  a  well-dclined  problem  exactly  like  that  posed  at  the  beginning 
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of  this  chapter;  therefore,  the  extended  modified  cut  metliod  can  again 
be  employed  as  follows:   place  the  largest  arc  in  A  -  A  -  F,  say  a  , 
into  the  set  A  and  solve  for  the  weakest  (n-| A| ) -disconnecting  set  in 
both  of  the  following  sets: 

V  =  {h  Ih   is  a  disconnecting  set  of  G(N,A-A)  ,  mn^^,  and,  a   e  H  } 
^    t   I  o  p    n 

and 
;V  ={h  |h   is  a  disconnecting  set  of  G(N,A-A)  ,  rnH=$,  and,  a   ^  H  }. 


The  General  Algoritlim 

The  general  algorithm  consists  of  3  routines.   Routine  A  is 
the  extended  modified  cut  method  with  the  set  T   of  infinite  capacitied 
arcs  and  a  set  A  of  arcs  to  be  removed  from  G(N,A)  (i.e.  ,  the  network 
G(N,A-A)  is  considered).   Routine  B  is  employed  when  Routine  A  ends 
in  a  success;  i.e. ,  no  gap  occurred  or  when  all  the  disconnecting  sets 
have  been  considered.   Routine  C  is  used  to  define  a  new  set  of  suIj- 
problems  if  a  gap  occurs. 

It  is  assumed  at  the  outset  that  A  and  T   are  both  empty. 
Initially  the  current  best  solution  value,  denoted  by  V  ,  is  set  to 
infinity  or  to  the  value  of  the  best  solution  known. 
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The  General  Algorithm  can  be  stated  as: 


ROUTINE  A 


STEP  0: 


If   n  -   A  =  0   let   u  =  oo,  otherwise   let   u  =  1. 


STEP   1: 


Let      c 


"(a)    =     ( 


;(a)      if      c(a)    <  u,    and      a   fi  T 
00        if      a   e  r 
u        if      c(a)    ^  u,    and      a  jt  V 


Step  2: 


STEP  3: 


Find  a  minimum  disconnecting  set  A   in  /O    (N,A-A).   Let 
U  =  {a  I  a  e  A   and   c(a)  =  u}. 


If   c"(A  )  =  00,  go  to  STEP  2,  ROUTINE  B. 


STEP  4: 


a.  If   |u|  >  n-JAl,  let   u  =  u+1   and  go  to  STEP  1. 

b.  If   |u|  =  h-IaI,  go  to  STEP  1,  ROUTINE  B. 
■c.   If   |u|  <  n-|A|,  go  to  ROUTINE  C. 
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ROUTINE  B 


STEP  1; 


STEP  2: 


If   C(A  -U)  <  V    let   A  =  A   and   V  =  C(A  -U) 
o  o 


If   A  =  $   terminate,  A   is  the  optimal  solution.   If 


A  5^  $   continue. 


STEP  3:' 


Let   A  =  A  -  [a  ]   where  p  ^  k  for  all   a   e  A 

(i.e. ,  a   is  the  smallest  capacity  arc  in  A). 

Let   r  =  r  U  {a  }  -  (r  n  {a   ,,a   ^, a  }) 

^  p-'         ••  p+1   p+2'      m-" 

(i.e. ,  add  a   to  F  and  remove  from  f  all  arcs  havine 
P  ^ 

indices  greater  than  p) .   Go  to  ROUTINE  A. 


ROUTINE  C 


STEP  1: 


STEP  2: 


If   A  -  A  U  r  =  $   go  to  STEP  2,  ROUTINE  B. 


Let  A  =  A  U  [a  ]  v.-here  p  S  k  for  all  a   e  A-A  U  F 
p  k 

(i.e. ,  a   is  the  largest  capacity  arc  not  in  A  or  F) 
Go  to  ROUTINE  A, 
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An  Improved  General  Algorithm 

In  this  section  results  which  improve  the  computational  effi- 
ciency of  the  general  algorithm  are  discussed  and  an  improved  general 
algorithm  is  developed. 

Lemma  4.1: 

If  a  cut  B   is  optimal  for  the  n-k  but  not  the  n  most  vital 

links  problem  of  G(N,A)  and  if  the  cut  H   is  optimal  for  the  n  most 

vital  links  problem,  then   C(H   -  H  )  >  C(B   -  B,  ) . 
'  n    k       n    k 

Proof: 

Note,  by  the  hypothesis  of  Lemma  4.1, 

C(B  -B,  )  S  C(H  -H  )  (2) 

k  k 

and 

C(B  -B  )  >  C(H  -H  ).  (3) 

n  n 

Combining  (2)  and  (3)  yields: 

C(B*-B*)  -C(B*-B*)  <C(H*-H*)  -  C(H-H  )  (4) 

k         n         k        n 

But  (4)  reduces  to   C(B  -B,  )  <  C(H  -H,  )  which  is  the  desired  result, 
n  k        n  k 

Q.  E.  P. 
Lemma  4.1  is  most  useful  for  the  case  where  k  =:  1 ;  that  is, 
consider  the  situation  where,  in  using  the  modified  cut  method,  a  gap 

occurs  at  the  n  but  not  at  the  n  -  1  most  vital  links  problem.   If  tlie 

*  th 

cut  B   is  found  to  be  optimal  for  this  latter  problem  and  the  n 

*  th 

largest  arc  of  B   is  a  ,  then  cuts  whose  n   largest  arc  is  no  larger 
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than  a    do  not  have  a  smaller  n-reduced  value  than  that  of  B   (the 

X 

optimal  cut  of  the  n-1  most  vital  links  problem). 

The  final  result  of  this  chapter  concerns  a  lower  bound  on 
the  optimal  solution  of  the  a   most  vital  links  problem  which  is  found 
as  a  by-product  of  the  modified  cut  metliod. 

Lemma  4.2: 

*  u  * 

If  B   is  a  minimum  cut  of  G  (N,A),  then  the  optimal  cut,  A  , 

to  the  a  most  vital  links  problem  satisfies: 

♦  *      u   * 

C(A  -A  )  >  C  (B  )  -  ucv  . 

Proof: 

Verification  of  this  resvilt  only  requires  the  observation  that 
if  B   contains  q'  u-ceiling  arcs,  then  B   is  the  optimal  cut  for  the  a 

most  vital  links  problem  (Theorem  3.2)  and  C(B  -B^,)  =  C  (B  )  -  ucv. 

*  u   * 
Thus,  the  optimal  solution  A   can  be  no  better  than  C  (B  )  -  ua. 

Q.  E.  P. 

Note  that  Lenma  4.2  also  holds  for  any  /O   (N,A)  where  one  is 
considering  all  cuts  but  those  which  contain  arcs  from  a  set  T    (this 
follows,  since  Theorem  4.2  and  Theorem  3.2  prove  parallel  results  for 
^    (N,A)  and  g"(N,A),  respectively). 

Although  given  gap  occurs,  the  bound  found  in  Lemma  4.2  is 
always  a  strict  lower  bound,  the  last  result  can  be  employed  in  two 
ways.   The  first  use  is  utilizing  the  bound  as  a  stopping  rule,  i.e. , 
if  the  current  best  solution  is  "close"  to  the  lower  bound,  terminate. 


54 


The  second  use  is  as  follows:   suppose  a  gap  occurs  at  some  subproblcm. 
If  the  lower  bound  on  the  optimal  of  this  subproblem  is  no  better  tlian 
a  solution  which  already  has  been  found,  then  the  subproblem  need  not 
be  considered  any  further,  even  if  it  ended  in  a  gap. 

Since  the  results  of  this  section  applj'  to  Routine  C,  an  imjiroved 
Routine  C  will  now  be  noted.   It  is  assumed  that  c(a  )  will  be  set  to 

X 

c(a  )  if  the  n  -  1  most  vital  links  are  not  found  when  solving  for  the 
m 

n  most  vital  links  in  G(N,A) ;  otherwise,  it  will  represent  the  capacity 

th 
of  the  n   largest  arc  of  the  optimal  cut  to  this  latter  problem. 

To  use  the  improved  algorithm  the  best  lower  bound  (BLB)  of  the  current 

subproblem  as  well  as  the  current  best  solution  and  the  current  best 

solution' s  value  (CBSV)  must  be  maintained. 

The  Improved  Routine  C  becomes: 

STEP  1: 

If   BLB  ^  CBSV,  go  to  STEP  2,  ROUTINE  B. 

STEP  2: 

Let  A  =  A  U  fa  ]  where  p  ^  k  for  all   a   e  A  -  A  U  T. 
^  p-^  k 

If   c(a  )  s  c(a  )   for  any   a   e  A,  go  to  STEP  2,  ROUTINE  B, 
q       X  q 


otherwise,  go  to  ROUTINE  A. 
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An  Example 


Recall  that  in  Chapter  3  it  was  establislied  that  a  gap  occurs 
in  solving  for  the  two  most  vital  links  of  Example  Netwoi'k  4.1, 
defined  by  Figure  4. 1.   The  Improved  General  Algorithm  will  now  be 
applied  to  this  network  to  resolve  this  gap  and  locate  the  two  most 


vital  links 


a  ,  100 


Figuz'e  4.1.   Example  Network  4.1. 


The  algorithm  starts  in  Routine  A  with  F  =  5,  A  =  $  (i.e. , 
the  original  network  G(N,A)).  As  noted  in  Chapter  3,  the  modified 

cut  method  ends  in  a  gap  with  the  following  minimal  cuts  of  G  (N,A) , 

1     2    r         1 
u  =  1,2,3,...,   being  located  A  =  A  =  i  ^^  >^9  >  ^-^  j  ^^^^ 

3  2      3 

a'  =  {a   ,a   ,a    a   }.   The  best  lower  bound  is  found  at  G   or  G 

A.  ^        JLo    J.  1    XO 

3   3 
and  has  value   C  (A  )  -  (2) (3)  =  8  -  6  =  2.   The  current  best  solution 

is  the  cut  [a   ,a   ,a    a   ]  with  a  2-reduced  cut  value  of  4. 

1.A  X*5    ii    XO 
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Since  a  gap  occurred,  Routine  C  is  employed  to  define  the 
network  depicted  in  Figure  4.2. 


33.100 


Figure  4.2.   ITie  First  Subnetwork  of 
Example  Network  4.1. 


Upon  definition  of  the  G(N,A-a  ) ,  Routine  A  is  again  employed 
to  locate  the  one  most  vital  link  of  this  network.   In  this  problem 

A"*  =  A^  rr  A   =  [Sg-^s^  ^^^      ^      =  ^^12 '^13 '^14 '^15^'   "^'^"^  ^  ^^^  occurs 
and  the  current  best  solution  value  remains  4.   The  best  lower  bound 
of  this  subprobleiu  has  a  value  of  4  and  is  found  at  iSf    (N,A).   Therefore, 
since  the  lower  bound  of  the  subproblem  is  not  better  than  the  best 
current  solution,  the  subproblem  does  not  have  to  be  considered  and 
Routine  B  is  called. 

Upon  applying  Routine  B,  the  network  sho\ra  in  Figure  4.3  is 
realized. 
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Figure  4.3.   The  Second  Subnetwork  of 
Example  Network  4.1. 


For  this  problem  the  cut  f ^^2 ' ^13 ' ^14 ' ^15^  ^^  *^®  minimum  cut 
for  every  ^    .       Thus,  at  J>    ,    the  two  most  vital  links  problem  is 
solved  for  the  above  network  and  Routine  B  is  called. 

Since   A  =  §,  the  algorithm  terminates  with  the  current  best 
solution  being  declared  the  optimal  solution  (i.e.,  any  pair  of  arcs 

from  the  set  [^^^2 '  ^13 '  ^14 ' '^'iS^  ^^^  ^^^   *^^°  ^°^^   vital  links  of  this 
network)  . 


CHAPTER   5 
COMPUTATIONAI.  RfJSULTS  AND  CONCLUSIONS 

The  modified  cut  method  of  Cliapter  3  and  the  improved  general 
algorithm  of  Chapter  4  were  coded  in  Fortran  IV  and  problems  were 
solved  using  an  IBM  360/65  computer.   These  algorithms  for  a  given 
network  require  approximately  the  same  amount  of  core  space  that  the 
maximum  flow  algorithm  of  Ford  and  Fulkerson  [12]  would  require   to 
solve  the  maximum  flow-minimum  cut  problem  for  the  same  network. 
For  the  example  networks  referred  to  in  this  chapter,  the  capacity  of 
the  arcs  ranged  from  1  to  50. 

Several  basic  network  configurations  were  used  to  test  the 
algorithms  and,  as  can  be  seen  in  Table  5.1,  the  test  problems  ranged 
in  size  from  small  (14  arcs,  6  nodes)  to  larger  networks  (250  arcs, 
50  nodes).   Also,  since  changing  the  capacity  function  of  a  network 
affects  tlie  solution  of  the  n  most  vital  links  problem,  several  of  the 
example  networks  have  the  same  structures  with  different  capacity 
functions. 


Actually  the  improved  general  algorithm  whicli  contains  the 
modified  cut  method  requires  four  )y\|- coraponent  vectors  more  than  the 
maximum  flow  algorithm. 

Naturally,  if  the  range  of  the  capacity  values  were  signif- 
icantly larger  than  50,  then,  for  a  given  large  problem,  one  might 
exj^ect  slightly  longei-  e.xecution  times  than  are  cited  in  this  chapter. 
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The  results  of  applying  the  modified  cut  method  to  the  example 
problems  are  presented  in  Table  5.1.   The  symbol  T]  which  has  been  used 
to  indicate  the  smallest  number  of  arcs  in  any  cut  between  tlie  source 
and  sink  nodes  in  G(N,A)  also  represents  the  ma.ximum  number  of  most 
vital  links  problems  that  could  be  of  interest  for  a  given  network, 
i.e.  ,  n  =  1,2,..  .  ,T\.      Those  values  of  n  for  which  a  gap  occurred  are 
noted  in  the  last  column  of  Table  5. 1  for  each  test  problem.   Observe 
that  for  the  sample  problems  which  were  solved, about  5  percent  of  the 
possible  most  vital  links  problems  ended  in  a  gap. 

In  Table  5.2  the  results  of  applying  the  improved  general 
algorithm  to  those  gaps  noted  in  Table  5.1  are  presented.   Since  the 
number  of  times  the  modified  cut  method  (Routine  A)  is  employed  is  a 
key  factor  in  determining  the  efficiency  of  the  improved  general  algo- 
rithm, this  statistic  is  included  in  column  6  of  the  Table  5.2. 


60 


TABLE   5.1 

RESULTS  OF  APPLYING  MODIFIED  CUT  METHOD 
TO  SAMPLE  NETWORKS 


Problem 
Number 

\a\ 

|n( 

^ 

Execution 

Time 
(minutes) 

Values  of  n 

for  Which 

Gaps  Occurred 

1 

28 

12 

4 

.03 

- 

2 

28 

12 

4 

.03 

- 

3 

28 

12 

4 

.03 

- 

4 

35 

16 

5 

.04 

- 

5 

35 

16 

5 

.03 

- 

6 

35 

16 

5 

.04 

- 

7 

14 

6 

4 

.03 

- 

8 

14 

6 

4 

.02 

- 

9 

14 

6 

4 

.03 

10 

52 

17 

8 

.07 

2 

11 

52 

17 

8 

.06 

- 

12 

52 

17 

8 

.  11 

4 

13 

41 

16 

8 

.06 

4,5 

14 

41 

16 

8 

.05 

- 

15 

41 

16 

8 

.07 

- 

16 

41 

12 

7 

.05 

5 

17 

41 

12 

7 

.03 

- 

18 

41 

12 

7 

.04 

- 

19 

48 

16 

7 

.04 

- 

20 

48 

16 

7 

.05 

- 

21 

48 

16 

7 

.08 
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TABLE    5.1    (Continued) 


Problem 

Number 

Execution 

Values  of  n 

|a| 

|n| 

T\ 

Time 
(minutes) 

for  Which 
Gaps  Occurred 

22 

26 

9 

7 

.05 

- 

•  23 

26 

9 

7 

.03 

- 

24 

26 

9 

7 

.05 

- 

25 

57 

20 

7 

.03 

- 

26  - 

57 

20 

7 

.05 

- 

27 

57 

20 

7 

.05 

- 

28 

35 

11 

8 

.04 

- 

2 

35 

11 

8 

.05 

- 

30 

35 

11 

8 

.05 

4,6 

31 

15 

7 

3 

.02 

1,2 

32 

250 

50 

25 

.30 

'.- 

33 

250 

50 

25 

,38 

6,18,19 

34 

250 

50 

25 

.35 

- 

35 

250 

50 

25 

.43 

6,7,8,9,12,13 

36 

250 

50 

25 

.36 

4,14 

37 

250 

50 

25 

.26 

- 

38 

250 

50 

25 

.50 

7,8,9,10,11 

39 

250 

50 

25 

.27 

- 

40 

250 

50 

25 

.29 

- 
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TABLE      5 .  2 


RESULTS   OF   APPLYING    IMPrxO\T::U    GENERAL  ALGORITHM 
TO   GAPS   THAT   OCCURRED    IN   SAMPLE   NETWORKS 


Problem 

Value  of  n 

Execution 

Times  Modified 

Number  in 

|a| 

|n| 

Solved  for 

Time 

Cut  Method 

Table  5. 1 

(minutes) 

Employed 

10 

52 

17 

2 

.48 

31 

12 

52 

17 

4 

.12 

3 

13 

41 

16 

4 

.09 

3 

13 

41 

16 

5 

.09 

3 

16 

41 

12 

5 

.07 

5 

30 

35 

10 

4 

.16 

11 

30 

35 

10 

6 

.08 

3 

31 

15 

7 

1 

.02 

3 

31 

15 

7 

2 

.02 

3 

33 

250 

50 

6 

2.72 

14 

33 

250 

50 

18 

1.10 

5 

33 

250 

50 

19 

1.06 

.5 

35 

250 

50 

6 

2.28 

7 

35 

250 

50 

7 

2.50 

7 

35 

250 

50 

8 

2.46 

7 

35 

250 

50 

9 

2.46 

8 

35 

250 

50 

12 

.92 

3 

35 

250 

50 

13 

.92 

3 

3(3 

250 

50 

4 

8.50 
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36 

250 

50 

14 

.69 

3 

38 

250 

50 

7 

12.15 

43 

38 

250 

50 

8 

12.20 

43 

38 

250 

50 

9 

12.20 

43 

38 

250 

50 

10 

12.56 

43 

38 

250 

50 

11 

12,04 

43 
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Conclusions 

To  find  the  n  most  vital  links  of  an  arbitrary  network  by 

enumeration  would  involve  solving  either  I      maximum  flow  problems 

\  n  / 

(one  for  each  possible  set  of  n  arcs  of  A)  or  locating  all  the  cuts 
between  the  source  and  sink  nodes.   As  noted  in  this  chapter  for  the 
problems  investigated  to  date,  a  gap  occurred  for  less  than  5  percent 
of  the  possible  n-values.   Thus,  95  percent  of  the  cases  required  one 
use  of  the  modified  cut  method  (which  is  effectively  the  same  as 
solving  a  maximum  flow  problem)  and  the  solution  found  included  a  fam- 
ily of  most  vital  links  of  which  the  n-value  of  interest  was  a  member. 
Furthermore,  even  in  situations  where  gaps  occurred,  on  the  average  it 
has  only  been  required  to  employ  the  modified  cut  method  sixteen  times 
per  gap  problem.   Tliis  is  approximately  equivalent  to  sixteen  maximum 
flow  problems. 

Next,  recall  that  two  bounds  have  been  identified  (Lemma  4.1 
and  Lenuna  4.2)  in  connection  with  the  improved  general  algorithm.   Wliile 
each  of  these  has  been  employed  in  particular  problems,  the  lower  bound 
on  the  value  of  the  optimal  solutions  (Lemma  4.2)  has  proved  to  be  much 
stronger  than  the  bound  noted  in  Lemma  4.1.   However,  since  both  bounds 
are  easy  to  find  and  maintain,  it  is  suggested  that  both  be  used  when 
applying  the  methods  developed  in  this  paper. 

In  conclusion,  while  the  above  observations  were  made  after 
considering  only  a  small  sample  of  problems,  the  partitioning  scheme 
of  the  improved  general  algorithm,  while  similar  to  implicit  enumer- 
ation (e.g.,  Balas  [1],  Geoffrion  [14]),  is  much  more  efficient  than 
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the  latter  ajjproach.   This  is  due  to  the  infrequency  of  gaps  occurring; 
the  bounds  which  truncate  subproblems ;  and   the  fact  that  each  time  an 
arc  is  placed  in  "  whole  classes  of  disconnecting  sets  are  eliminated 
from  fui'ther  consideration. 


APPENDIX  A 
FORMULATIONS  OF  THE  n  MOST  VITAL  LINKS  PROBLEM 

A  Linear  Formulation 

In  Chapter  3,  it  was  demonstrated  that  the  n  most  vital  links 
problem  could  be  posed  as  a  cut  problem  where  the  optimal  cut  A* 
satisfies 

C(A  -A  )  =  min   C(B-B  )      for  all  cuts  B  (A. 1) 

n  n 

and  A   contains  the  n  most  vital  links.   This  cut  problem  will  now 
n 

be  formulated  as  a  linear  zero-one  integer  programming  problem.   First 
observe  that,  instead  of  only  considering  cuts  in  equation  (A. 1) ,  the 
B's  could  represent  all  disconnecting  sets  separating  the  source 
and  sink.   This  is  valid,  since  the  weakest  disconnecting  set  will 
always  be  a  cut. 

The  minimum  cut  problem  as  formulated  by  Ford  and  Fulkerson  [11] 
may  be  posed  as  a  zero-one  integer  programming  problem,  i.e. , 

Min   E   c(i,j)  Y(i.J) 
(i.j) 

such  that  -rr(s)  +  ■n(t)  2:  1 
I 

TT(i)  -  TT(j)  +  Y(i.J)  -  0      (i,J)  e  A 

TT(i),Y(i'J)^0.1    (i,j)eA  and  (i)  e  N 
65 
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Furthermoi-e,  a  solution  is  feasible  for  I,  if  and  only  if,  the  set 
of  arcs  for  which  7(1. J)  =  ^  forms  a  disconnecting  set  of  the  given 
network. 

Now  define  a  now  set  of  zero-one  variables  u  such  that  for  each 
arc  (i.j)  there  is  a  variable   u(i,j).   These  u(i,j)'s  are  restricted 
to  satisfy   u(i,j)  ^  y(i,j)-      Thus,  for  any  solution  only   u(i,j)'s 
corresponding  to  arcs  in  the  identified  disconnecting  set  may  be  non- 
zero  Adding  the  constraint   I   u(i,j)  ^  n  insures  that  no  more  than 
n  of  the  u(i,j)  variables  will  be  non-zero  for  any  feasible  solution  of 
the  problem  to  be  defined.   The  n  most  vital  links  problem  can  then  be 
formulated  as: 

Min   E   c(i,j)  Y(i.J)  "    ^   c(i,j)  u(i,j) 

(i,j)  (i.j) 

such  that 

-TT(s)  +  TT(t)  s  1      (i,j)  e  A 

^^  TT(i)  -  TT(j)  +  Y(i'J)  ^  ^  (^'J)  ^  ^ 

,.2  u(i,j)  S  n       (i,j)  e  A 
(i.j) 

u(i,j),  ydj),  TT(i)  =  0,1       (i,j)  e  A   and   (i)  e  N 
u(i,j)  -  Y(i.J)  ^  0        (^-J)  ^  ^ 
Note  that  only  disconnecting  sets,  B,  are  feasible  for  II  and 

for  each  such  B  the  optimal  u(i,j)  values  are  one  if  (i,j)  e  B^  and 

zero  otherwise.   Thus,  the  optimal  solution  to  II  is  that  disconnecting 

set  A  which  satisfies  equation  (A. 1) . 
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By  making-  the  change  of  variable  y  (i»j)  =  Y^-'-'J^  ~  '^(^iJ)' 

problem  II  becomes: 

Min   E   c(i,j)  y    (i,j) 
(ij) 

such  that 

-Tt(s)  +  Tr(t)  ^  1 
III 

Tr(i)  -  rr(j)  +  y'^i-j)  +  "(i.J)  ^  0     (i.j)  ^  A 

E  u(i,j)   s  n     (ij)  e  A 
•  -  (ij) 

u(i,j),  y'i±,3),    TT(i)  =  0,1    (i,j)  e  A  and  (i)  e  N 

A  Quadratic  Formulation 

In  this  section  a  zero-one  quadratic  formulation  with  special 
structure  is  developed  for  solving  the  n  most  vital  links  problem. 
Basically,  the  formulation  is  Formulation  II  presented  in  the  last 
section  with  the  constraints   y^^'J)  ~  u(i,j)  ^  0   removed  from  the 
constraint  set.   Tliese  conditions,  i.e.,  u(i,j)  ^  y(i,j),  are  guar- 
anteed at  optimal ity  in  the  new  formulation  by  introducing  a  quadratic 
part  to  the  objective  fimction  which  is  zero  if   u(i,j)  ^  y(i,j)  and 
CO  if  the  constraint  is  violated.   The  objective  function  will  be  shown 
to  be  concave  and  the  constraint  coefficient  matrix  unimodular.   ITiis 
structure  will  allow  the  optimal  solution  to  this  new  formulation  to 
be  found  by  ranking  the  extreme  points  of  a  related  continuous  quad- 
ratic programming  problem. 
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Consider  the  following  zero-one  quadratic  programming  problem, 
where  M  is  some  very  large  positive  number. 


Min   E   c(i,j)  y^i.J)  "   ^   c(i,j)  u(i,j) 
(i.j)  (i,j) 

2  2l 

+   E  m{i  -^(u(i  J)  -Y(i,j)  -  l)  -|(u(i,j)  +Y(i-J)  -  l)  '' 
(i  ,j) 

such  that 

-n(s)  +  TT(t)  S:  1 

TT(i")  -  TT(j)  +  Y(i>J)  ^  0      (i'J>  ^  A 


IV 


E   u(i,j)  =2  n     (i,j)  e  A 
(i  J) 

u(i,j),  Y(i-J)'  ^(i)  =  0-1    (i'J^  ^  ^   ^"^^   (i)  ^  N 


Let  f(^u(i,j)  ,Y(i,j))  =   {l-^^uCi.j)  -yCiJ)  -  l)^ 

2 

-|(u(i,j)+Y(i.j)-l)  } 


2 

'  M 


Note  that   f  I  u(i  ,  j)  .yd  ,  j  )  )  =  0  whenever   yCiJ)  >  u(i,j)  in  IV  and 
ffu(i,j) ,y(i,j)  )  =  M  if   y(i,j)  =  0   and  u(i,j)  =  1.   Thus,  for  M 
sufficiently  large,  the  optimal  solutions  to  Formulation  II  of  the  last 
section  and  IV,  given  above,  will  be  the  same  and  IV  will  solve  the 
n  most  vital  links  problem. 

Next,     E   f (u(i , j) ,y(i , j)  j  will  be  shown  to  be  a  concave 

(ij)  \  X 

function  of  u    and  \^   .      Note  that   f  i  u(i  ,  j  )  ,y(i  ,  j)  1  is  a  sum  of 
a  constant  and  two  quadratic  forms,  each  of  which  has  a  value  of  less 
than  or  equal  to  zero  for  all  values  of   y(i,j)  and  u(i,j);  thus,  each 
of  the  forms  are  concave.   Furthermore,  the  sums  of  concave  functions 
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are  concave  so  f ( u(i , j ) jyCi , j)  ]  is  concave  as  is  the  function 

E   f(Y(i,j),u(i,j) 
(ij) 

Consider  Formulation  V  which  is  really  a  continuous  version 
of  Formulation  IV. 


Min   E   c(i,j)  yii.j)    -        2   c(i,j)  yCi.J)  +   ^  f ( u(i , j ) .yCi , j) 
(i,J)  (iJ)  (iJ) 

such  that 

-tt(s)  +  rr(t)  s  1 
V 

TT(i)  -  TT(j)  +  Y(i.J)  ^  0     (i'j)  ^  A 

E  u(i,j)  S  n     (i,j)  e  A 
(iJ) 

0  ^  Y(i.j)  ^  1      (i.j)  e  A 

0  ^  u(i,j)  s  1     (i,j)  e  A 

0  S  TT(i)  2  1       (i)  e  N 

Hadley  in  [16]  notes  that  when  minimizing  a  concave  function 
over  a  convex  set,  there  is  an  optimal  to  the  problem  which  occurs 
at  an  extreme  point  of  the  convex  set.   Tlie  convex  constraint  set  of 
V  will  now  be  shown  to  be  unimodular.   Thus,  since,  if  this  is  the 
case,  each  extreme  point  solution  will  be  zero-one,  the  optimal 
zero-one  solution  to  IV  can  be  found  by  ranking  the  extreme  points 
of  the  constraint  set  of  V.   Murty  [23],  and  Cabot  and  Francis  [8] 
demonstrate  how  to  rank  the  extreme  points  of  quadratic  programming 
problems. 
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It  only  remains  to  verify  that  the  coefficient  matrix  of  V 
has  the  unimodularity  property.   The  first  two  sets  of  constraints  of  V 
are  just  the  constraints  of  the  minimum  cut  problem  and  it  is  a  known 
fact  that  the  coefficient  matrix  associated  with  the  constraint  set  of 
the  minimum  cut  problem  is  unimodular.   Tlius,  by  letting  the  coefficient 
matrix  associated  with  these  constraints  be  denoted  by  R  the  problem 
reduces  to  showing  that  if  R  is  a  unimodular  matrix,  then  R   is  also 
luiimodular. -  Where  R   is  defined  in  partitioned  form  as: 


R   = 


R 

1 

0 

0 

1 

1 

0 

1 

I 

where  0,  1  are  row 
vectors  and  a   is  equal  to 

2|aJ  +  |n|  for   G(N,A) 


The  matrix  R  is  unimodular  if  and  only  if  every  minor  determinate  is 

equal  to  0,  +1,  or  -1.   To  show  this  result  holds  for  R  ,  the  matrices 

R,  ,  R„,  R^  will  be  defined  where  R,  =  [R:0],  R„  r.  [0:1]  and  R^  =  [0:1  ]. 
Iz3  1  2—  —        JQ' 

Let  A  be  an  arbitrary  minor  of  R  .   If  A  fl  R  ;t'  $,  then  there  is 

some  row,  r,  of  A  which  contains  elements  of  R  .   By  expanding  A  by  this 

row  a  matrix,  B,  satisfying   B  s;  A  and  |A|  ::r  ±  j  B|  will  be  found. 

Also,  B  will  not  contain  the  row  r.   Likewise,  it  is  possible  to  exjiand 

B  by  any  row  whicli  contains  elements  of  R   and  this  expansion  process 

can  be  continued  until  for  some  matrix,  C,  such  that  \a\    =  ±  |c(,  and 


c  n   R  =  §. 
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Next,  consider  the  case  where  a  minor,  C,  of  R  has  C  H  R„  =  $. 
If   C  n  R  =  $,  then   C  c  R      and  thus   | c|  =  0 ,  +1  or  -1 ,  since  R   is 
unimodular.   If   C  n  R  ?^  $ ,  then  there  are  colunuis  of  C  which  contain 
elements  of  R^ .   By  a  sequential  expansion  about  this  set  of  columns, 
it  will  be  possible  to  equate   |c|  =  ±  |d[  where   D  C  C   and   D  0  R2  =  ^ 
But,  Dec  implies  that  D  0  R^  =  $  which  in  turn  implies  that   D  c  R 
Thus,  |c|  =  ±  (d|  which  has  a  value  of  0,  +1,  or  -1. 

The'  above  discussion  leads  to  the  following  method  for  showing 

that  if  A  is  a  minor  of  r',  then  |a|  =  0,  +1,  or  -1.   First  A  is 

expanded  by  any  row  containing  elements  of  R  and  the  method  described 

above  is  used  to  get  a  matrix  B  where:   B  C  A,  \a\^   ±    \b\  ,    and  B  D  R„ =  ^ 

Next,  the  columns  of  B  are  inspected  to  find  if  B  H  R  =  $.   If  this  is 

the  case,  B  c:  R  ,  and  the  desired  result  has  been  found.   If   R  (1  R  7^  ? 
1  ^ 

then  the  above  expansion  scheme  will  locate  a  new  matrix,  C,  Such  that 
(a)  :=  ±  |b|  =  ±  |c|  and   C  gR  .   So  in  either  case   |b|,  and  thus  |a|, 
has  a  value  of  0,  +1,  or  -1. 
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